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Chapter  1 :  Introduction 

Chapter  1 ■ t  Comparison  between  proposed  and  achieved  tasks 

The  aims  of  our  study  as  outlined  in  our  original  proposal  were  the  following: 

(a)  Carry  out  an  initial  study  of  suitably  chosen  LANDSAT  data  over  both 
land  and  sea.' 

(b)  Develop  the  first  step  in  the  evolution  of  semi-automatic  computer 
algorithms  to  obtain  cloud  morphology  data  on  a  global  and  long-term 
bas  i  s 

(c)  Provide  a  basis  for  future  coordinated  satellite,  air  and  ground-based 
experiments  on  cloud  field  properties. 

(d)  Test  and  improve  the  modelling  of  cloud  field  optical  properties  for 
various  purposes.' 

(e)  Use  the  data  for  testing  models  of  atmospheric  dynamics  on  various 

scales  and  to  develop  methods  for  using  cloud  morphology  data  as  diagnostics 
for  large  scale  atmospheric  characteristics. 

These  aims  have  been  carried  out  for  the  most  part. ^ Sui table  LANDSAT  2,3 
and  1*  imagery  was  chosen  over  the  Atlantic,  Pacific  and  Indian  Oceans.  Hundreds 
of  images  were  examined  by  eye  and  computer  catalogue  at  the  Goddard  Space  Flight 
Center. 

Thirty  LANDSAT  2  or  3  full  scenes  and  41  LANDSAT  4  quarter  scenes  were  selected 


of  which  only  the  LANDSAT  2  and  3  scenes  have  been  analysed.  The  LANDSAT  4  images 
which  include  thermal  IR  data  are  in  the  process  of  being  analysed  and  will  be 
the  subject  of  a  separate  report.  This  will  also  include  the  calibration  and  use 
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of  a  new  method,  based  on  the  solar  spectral  range,  of  cloud  top  height  deter¬ 
mination,  developed  by  us  and  described  in  previous  reports. 

We  have  developed  a  semi -automated  algorithm  for  the  production  of  cloud 
field  maps  over  sea.  Because  of  the  unexpected  degree  of  inhomogeneity  of  the 
ocean  surface,  we  studied  a  variety  of  schemes  for  cloud-background  discrimination 
none  of  which  is  really  optimal.  This  is  also  the  reason  the  analysis  of  imagery 
over  land  has  been  deferred  for  a  later  study. 

The  method  has  been  applied  and  is  suitable  for  man-computer  interactive 
analysis  of  cloud  fields  after  suitable  preparation  of  the  analyst.  The  method 
provides  hi stogrammes  of  single  cloud  parameters  and  their  various  moments  and 
correlation.  Some  ideas  have  evolved  for  the  analysis  of  2D-cloud  patterns  which 
have  been  submitted  to  the  USAF  in  a  "letter  of  intent"  recently. 

The  methods  of  analysis  and  the  data-base  may  be  used  as  a  basis  for  the 
planning  of  possible  future  dedicated  satellite-based  cloud  field  morphology 
missions.  These  should  include  surface  and  air-based  field  morphology  data  at 
all  t imes  of  the  day . 

We  have  developed  an  analytical  model  of  the  reflectivity  of  bar  clouds 
of  given  aspect  ratios,  inter-cloud  distances  and  optical  properties.  Cloud- 
cloud  scatter  and  shading  have  been  explicitly  taken  into  account.  Future 
development  will  include  a  3-D  box  cloud  array  and  the  use  of  our  data  base  for 


cal i brat  ion . 
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Finally,  the  data  base  has  been  applied  to  a  current  model  of  cloud  convection 
and  mean  convective  mass  fluxes  have  been  estimated  which  are  in  reasonable 


agreement  with  theoretical  and  in-si tu  estimates. 

1.2:  Presentations  at  scientific  meetings 

The  materia)  summer i zed  in  this  report  has  been  presented  in  part  at 


several  meetings. 


J.H.  Joseph:  "The  Morphology  of  Broken  Cloud  Fields  from  LANDSAT". 
Symp.  R1  -  Interdependence  of  clouds.  Trace  Substances,  Radiation  and 
Climate  IAMAP  -  IAGA  Conf.  Hamburg,  F.R.G.  August  1983. 


2. 


J.H.  Joseph:  "The  Morphology  of  Fair  Weather  Cumulus  Clouds  as  Remotely 
Sensed  from  Satellites  and  some  Applications"  Paper  4.3.**.  Symp.  on  Space 
Observations  for  Space  Studies,  XXV  COSPAR,  Graz,  Austria,  25-29  June  1984. 


3.  R.F.  Cahalan:  "Climatological  Statistics  of  Cloudiness"  -  Invited  Paper, 
Proc.  5th  Conf.  on  Atmospheric  Radiation,  Baltimore,  MD,  Oct.  1983. 


J.H.  Joseph:  "The  Role  of  Cloud  Field  Morphology  in  Weather  and  Climate 
Studies  -  The  Charney  Workshop  on  Weather  Modification,  Caesarea,  Oct.  11-13, 
1982. 


The  material  will  be  discussed  in  detail  at  the  IBM  Workshop  on  "Fractals  and 
Dynamics"  at  Oberlach,  Austria,  June  30  -  July  7,  1986. 


I .  3  Publ  i  caTi  ons 

The  following  publications  have  resulted  from  our  study: 

J. H.  Joseph  (1985):  The  morphology  of  fairweather  Cu  clouds  as  remotely  sensed 
for  satellites  and  some  applications,  Adv.  Space.  Res.  5,  213-216. 

J.H.  Joseph  (1985):  The  role  of  cloud  field  morphology  in  weather  and  climate 
studies,  Is r.  J.  Earth  Sci.  34,  96-101. 

-N 

J.H.  Joseph  and  R.F.  Cahalan:  On  the  morphology  of  broken  Cumulus  cloud 
fields  (submitted)  1 986 . 


The  various  chapters  in  this  report  essentially  represent  the  papers  that 

will  be  submitted  for  publication  in  the  near  future. 

1 . 4  Future  Applications 

The  application  of  our  results  to  the  calculation  of  transmission  of 
radiation  or  information  through  cloudy  layers  for  various  purposes  -  e.g.  correction 
of  remotely  determined  sea  surface  temperatures,  availability  of  solar  radiation 
or  the  probabi 1 i ty.  of  lines  of  sight  will  be  discussed  this  fall  at  the  NASA  - 
GSFC,  Greenbelt,  MO.,  which  I  plan  to  visit  for  an  extended  stay. 


y  '■y-M'ivf  y-.y.vJtvX  ’ 
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1 .5  Work  in  Process 


This  report  does  not  Include  the  totality  of  tasks  outlined  in  our 
proposal  because  the  USAF  funding  was  discontinued  prematurely. 

We  are  in  the  process  of  analysing  a  large  amount  of  LANDSAT  k  TM  data 
and  have  an  option  to  acquire  SPOT  data.  The  rate  of  analysis  is  slower  now 
because  of  local  improvised  funding. 

The  analysis  planned  is  similar  and  we  expect  to  achieve  the  following 

alms : 

(a)  Improve  our  resolution  of  cloud  sizes,  etc,  down  to  30  meters; 

(b)  Add  an  additional  source  of  information,  namely  the  thermal  IR 
emissivi ties; 

(c)  Determine  the  distribution  of  the  heights  of  cloud  tops  from  the 
emission  temperatures  using  spatial  coherence  techniques.  These 
heights  for  the  first  time  will  be  compared  to  those  calculated  with 
our  model  for  height  determination  in  the  solar  range,  using  data 
from  the  same  satellite  platform; 

(d)  Get  a  time  series  on  a  monthly  or  yearly  basis  of  possible  variations 
in  cloud  properties  primarily  over  the  ocean  East  and  West  of  Florida. 
This  will  make  it  possible  to  correlate  the  latter  with  year  to  year 
changes  in  meteorological  conditions. 


13  - 


Chapter  2:  Processing  of  LANDSAT  DATA 
1 .  Introduction 

The  type  of  LANDSAT  Imagery  analysed  in  this  study  is  broken  Cumulus 
cloud  fields  over  ocean. 

The  reasons  are  several  and  are  summarised  here  in  short  for  complete¬ 
ness.  Firstly,  about  50%  of  all  clouds  are  of  this  type  and  exist  almost 
continuously  over  great  parts  of  the  globe.  Cumulus  clouds  are  therefore 
an  important  part  of  the  global  cloudiness  and  are  thus  worthy  of  investiga¬ 
tion  in  their  own  right. 

Secondly,  as  will  be  described  in  this  chapter,  the  discrimination 
of  clouds  against  their  background  is  one  of  the  biggest  problems  of  the 
data  analysis.  Bright  clouds  over  a  dark  sea  were  therefore  chosen  as  a 
first  topic  of  invest igat ion ,  later  augmented  by  a  study  of  the  effects  of  a 
quite  variable  ocean  background  on  the  results  of  the  analysis. 

Four  main  classes  of  problems  will  be  addressed  in  this  chapter: 

(1)  Discrimination  between  cloud  and  background 

(2)  Recognition  of  a  connected  cluster  of  pixels  as  a  cloud 

(3)  Definition  of  the  properties  of  a  single  cloud 

(a)  area; 

(b)  perimeter; 

(c)  effective  radius; 

(d)  location; 

(e)  center  of  gravity; 

(f)  nearest  neighbour; 

(g)  nearest  neighbour  of  comparable  or  larger  area; 

(h)  holes  in  cloud; 

(i)  cloud  top  altitude. 
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(4)  Properties  of  the  cloud  field,  assuming  clouds  to  be  independent 
of  one  another.  Both  the  experimental  data  and  their  fit  to 
chosen  functions  are  available. 

(a)  Distribution  function  of  equivalent  radii,  rft; 

(b)  Distribution  function  of  areas,  a; 

(c)  Distribution  function  of  nearest  neighbours.  D; 

(d)  Distribution  function  of  perimeters,  P; 

(e)  The  two-dimensional  histogram  of  (log  P,  log  a)  and  the 
fractal  dimension  of  cloud  fields; 

(f)  The  relationship  of  actual  cloud  area  to  area  of  influence 
of  a  cloud; 

(g)  Total  cloudiness; 

(h)  Total  number  of  clouds; 

(i)  Area  of  analysed  image 


2. 


Discrimination  between  cloud  and  background 


The  discrimination  between  cloud  and  ocean  background  is  not  a  simple  task. 

Even  though  brightness  of  the  ocean  surface  is  often  assumed  to  be  relatively 
homogenuous,  this  is  not  the  case  in  actuality.  The  brightness  of  ocean  surface 
changes  because  of  variations  in  depth  and  the  presence  of  waves,  haze  or  fog. 

In  order  to  determine  the  location  of  the  cut  between  cloud  and  background, 
we  applied  or  developed  several  methods,  no  single  one  of  which  turned  out  to  be 
optimal  by  itself. 

2 . 1  The  incremental  method 

This  method  starts  by  assigning  a  lower  limit  to  cloud  brightness  (gray 
level)that  is  much  too  high.  One  then  calculates: 

(a)  the  total  number  of  clouds; 

(b)  the  number  of  clouds  less  than  4  pixels  in  sizej 

(c)  the  total  relative  cloud  cover 

Next,  lower  and  lower  thresholds  are  assigned  and  the  above  three  properties  are 
calculated  for  each  case.  The  process  terminates  when  one  or  more  of  the 
following  phenomena  occurs 

(a)  the  number  of  small  clouds  increases  steeply.  This  means  that 
background  noise  is  added  to  the  actual  cloud  pixels; 

(b)  the  number  of  clouds  of  area  larger  than  4  pixels  decreases 
significantly  or  again  random  noise  connects  large  clouds; 


(c)  the  percentage  of  cloudiness  increases  s ign i f i can t ly . 


-16- 

Finally,  we  assume  that  we  have  converged  on  an  appropriate  threshold  for  a 
given  Image  when  two  different  spectral  bands  give  the  same  number  of  clouds  and 
the  same  cloudiness 

The  problem  with  this  method  is  that  it  is  difficult  to  determine  a  given 
fixed  rate  of  increase  of  the  small  clouds  or  rate  of  decrease  in  the  number  of 
large  clouds  chat  fits  all  available  Images. 

Experience  has  shown  that  these  criteria  change  from  image  to  image  and  even 
from  spectral  band  to  spectral  band  in  a  single  image. 

Figure  2.1  and  Table  2.1  show  a  typical  example  of  the  application  of  this 
method  in  case  7025. 

Figure  2.1  shows  the  increase  with  decreasing  values  of  the  chosen  threshold 
of  the  grey  level  of  the  number  of  "all  clouds"  (Figure  2.1(a)  -  MSS4,  Figure  2.1(f) 

MSS7), of  the  "big  clouds  (Figure  2.1(b)  -  MSS4,  Figure  2.1(g)  -  MSS7),  of  the  "small 

clouds"  (Figure  2.1(c)  -  MSS4,  Figure  2.1(h)  -  MSS 7),  of  the  amount  of  cloudiness 
due  to  the  "big  clouds"  (Figure  2.1  (d)  -  MSS4,  Figure  2.1(i)  -  MSS 7)  and  of  the 

amount  of  cloudiness  due  to  "all  clouds"  (Figure  2.1(e)  -  MSS4,  Figure  2 . 1 ( j )  -  MSS7) 

All  parameters  increase  slowly  with  decreasing  values  of  the  threshold.  Because  of 
the  slow  variations,  we  show  the  same  data  once  more  in  Table  2.1.  On  analysis  of 
the  figures  and  the  Table  the  reason  for  the  actual  choice  of  optimal  thresholds  for 
the  analysis  of  69  for  MSS4  and  of  2  for  MSS7  becomes  clear. 

2 . 2  The  composite  histogram  method 

When  developing  this  method  we  tried  to  set  up  an  automatic  way  to  exclude 
the  background  by  using  composites  of  appropriate  sample  histograms. 

Let  us  assume  that  the  general  histogram,  n(x),  includes  both  cloud 
and  background  pixels  at  any  given  brightness  or  gray  level  x.  Then  the 
probability  for  a  pixel  to  have  the  gray  level  x  is 

P  (x)  =  n  (x)/N 

where  N  is  the  total  number  of  pixels  in  the  histogram.  Assuming  we  have 
a  mixture  of  cloud  and  background  pixels  independent  of  each  other,  then 
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m,  =  N  /N 

D  D 


pb(x)  ■  nb(x)/Nfa 

P  (x)  -  n  (x)/N 
c  c  c 


and  "b"  and  “c"  denote  "background"  and  "cloud"  respectively.  In  addition, 


m.  +  m  =1 
b  c 


Next,  let  us  assume  that  for  brightness  lower  or  eaual  to  the  maximal 
brightness  of  background,  x  ,  there  are  no  cloud  pixels,  or  for 


x  <  x  ,  :  P  (x)  =  0 
—  maxb  c 


This  assumption  will  often  but  not  generally  be  true  for  clouds  over  the 
ocean,  but  not  necessarily  so  over  land.  We  thus  have  in  the  region 

(  X  —  Xmaxb  ) 

P(x)  =  pbM  mb 


In  particular  at  x  ■  x 


maxb 


m. 


P(xmaxb,/  Pb(xm»b) 


n(xma»b>/W 

"b(xmaxb,/Nb 


Therefore  the  probability  distribution  function  for  cloud  brightness  as  a 
function  of  x  is  given  by 

,  ,  P(x)  -  P,(x)  n, 

C  «  - : - 

I  -  mb 

The  parameter  m^  is  assumed  to  be  independent  of  x  and  determinable  from 
suitably  chosen  samples  from  a  given  image. 


In  practice,  this  method  is  applied  in  the  following  manner.  Using 
a  film  positive  in  a  given  band,  n  samples  of  cloud-free  areas  and  m  samples 
of  partly  cloudy  areas  are  chosen  by  eye. 


A  probability  distribution  function  for  the  background  is  constructed 
from  the  digital  data 


n  n 

P.  (x)  -5  HI .  *  Z  H.  ./  Z 

'  '  j.i  '•Jj.. 


A. 
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where  A.  is  a  sample  area  in  pixels  and  H.  . 

J  i  .J 

background  sample. 


is  the  histogram  of  a 
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v ,v.y.y "■ 


Similarly,  the  probability  distribution  function  for  the  mixed  samples 


is  constructed. 


P£(x)  i  HI. 


m  m 

I  H !  ./  Z 

j-1  *’J  j-1 


A. 

J 


where 


A*. 

J 


and  h!  . 
i  ,J 


are  the  area  and  histogram  of  a  mixed  sample. 


The  probability  distribution  function  for  the  clouds  alone  is  then 
given  by 

P3(x)  =  H3.  =  (H2.  -  H 1  x  R)/  (1  -  R) 

whe  re 

R  =  m.  -  H2.  /HI. 

b  i  i 

max  max 


The  index  i  is  taken  as  that  i  where  HI.  reaches  its  maximal  value  and 
max  i 

H3.  is  put  equal  to  zero  wherever  it  is  negative. 

Because  of  the  possible  uncertainties  in  the  application  of  this  method 
we  determine  an  intermediate  range  of  gray  levels,  in  which  both  clear  and 
cloud-filled  pixels  may  exist.  We  define  this  range  to  lie  between  the  lower 
limit  for  the  existence  of  cloud  pixels  and  the  lower  limit  for  "doubtful" 
pixels.  The  pixels  in  this  range  are  included  with  the  clouds  when  they  are 
next  to  clouds  and  disregarded  when  they  occur  alone. 


E 


The  i ntersect ron  between  the  cloud  and  background  distribution  functions 

is  defined  as  the  lower  limit  for  cloud  pixels  and  denoted  ICLOUD.  It  is 

determined  from  the  data  in  the  following  manner.  The  difference  function 

0.  -  H3,  -  HI. 

is  scanned  from  high  to  low  grey  level.  When  the  first  negative  value  is  obtained, 
we  calculate 

min  (0. ,  D.+]) 

and  ICLOUD  is  set  as  the  gray  level  with  the  lowest  difference,  D. • 

Next,  in  order  to  define  the  doubtful  range  -  IQUES  -  we  again  scan  from  high 

to  low  starting  from  the  maximum  of  H3. ,  looking  for 

H3.  <  -  max  (H3) 
i  —  e 

The  gray  level  that  fulfils  the  above  condition  is  defined  -  tentatively  - 
as  IQUES  1.  If 


IQUES  J  ?■  ICLOUD, 

we  look  for 

H3.  £  — ^  max  (H3) , 
e 

which  is  then  defined  as  IQUES  1. 


Finally,  we  re-scan  the  difference  function  from 


I MIN  -  min  (I CLOUD,  max(H3)) 

to  find  the  first  grey  level  for  which  the  function  becomes  zero  or  negative.  This 
point  is  denoted  IQUES  2  and  the  lower  limit  of  the  doubtful  range,  IQUES,  is  taken  as 

IQUESS  -  max  (IQUES  1,  IQUES  2) 

Figure  2.2  for  case  7025  shows  H3,  H2  and  HI  as  well  as  the  various  grey  levels 
described  above. 

For  the  case  of  MSS7  this  method  gives  a  threshold  value  similar  to  the  one 
obtained  by  means  of  the  previous  method  2  and  moreover  one  can  see  a  separation 
between  the  dotted  histogram  for  the  clouds  alone  and  the  full-line  curve  representing 
the  ocean  background.  This  is  not  the  case  however  for  MSS4  where  the  two  histograms  an 
poorly  separated  and  the  previous  method  indicated  a  working  threshold  of  69  as  opposed 
to  the  value  of  55  arrived  at  here. 

The  last  operation  carried  out  is  an  attempt  to  assess  the  effect  of  the 
particular  choice  of  samples  made  in  the  above  process. 

We  evaluate  a  standard  deviation  for  each  gray  level: 

S.  =  /  £  .  (H.  .  -  M.)Z/n 

i  J=1  i,J  i 

where  M.  is  the  average  value 

n 

M.  =  Z  H.  ./n 
'  j=]  '*J 

Two  new  probability  distribution  functions  are  constructed  for  the  background 

n 

HIP.  =  (Z  H.  .  +  S . ) / AP 


Where  AP  and  AM  are  the  histogram  areas  to  which  the  area  due  to  the 


standard  deviation  is  added  or  subtracted  respectively. 

One  then,  once  more,  goes  through  the  calculation  of  the  probability 
distribution  function  for  clouds,  namely  H3Pi  and  H3M  ,  for  the  two  perturbed 
background  functions,  HIP^,  HlM^.  Then  once  more,  one  evaluates  the  areas  of 
uncertainty  through  IQESP ,  ICLOUDP,  IQUE3M  and  ICLOUDM. 

In  Figure  2.3  we  show  the  function  HIP,  H3P,  HIM,  H3M  showing  the  possible 
error  due  to  the  sampling  to  be  unimportant. 

This  method  gives  good  results  for  low  gain  data  and  does  not  work  for 
high  gain  data. 


2 . 3  Spatial  coherence 

The  spatial  coherence  technique  is  based  on  an  analysis  of  the  dependence 
of  the  variance  of  the  intensity  on  the  latter  in  a  given  image.  This  is  done 
by  choosing  an  area  and  dividing  it  into  rectangles  of  nm  pixels.  In  each 
rectangle,  one  calculates  the  average  gray  level  and  its  standard  deviation. 

These  latter  are  plotted  versus  the  average  intensity  and  then  one  looks  for 
regions  of  low  standard  deviation  or  high  spatial  coherence  of  the  pixels,  ("feet") 
Normally,  one  should  get  two  regions  of  low  standard  deviation  -  one,  for  low 
intensity,  representative  of  the  ocean  background,  the  other,  for  high 
intensities,  represents  fully  cloud-covered  pixels  in  regions  of  thick  clouds. 
Between  these  "feet"  one  finds  an  arch  of  points  representing 
pixels  with  a  cloud-ocean  mix  of  intensities  or  covered  by  thinner  clouds. 
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In  our  case,  however,  the  method  does  not  work  very  well.  We  do  get  a 
narrow  cloud  "foot"  or  "double  foot".  The  foot  representing  the  ocean  is  very 
wide  and  is  not  suitable  for  discriminating  between  clouds  and  ocean.  In  other 
words,  the  ocean  surface  is  not  very  homogenuous. 

When  analysing  the  thermal  IR  data  of  LANDSAT  4  and  5,  the  method  does  work 
and  will  be  useful. 


2 . 4  Multispectral  histograms 

Using  MSS4  and  MSS 7,  we  constructed  one  and  two-dimensional  histograms  in 


both  digital  and  graphical  form. 


Histograms  were  constructed  of  the  whole  image,  of  cloud  filled  and  of 
clear  background  areas.  The  various  histograms-  and  the  images  were  then  analysed 
visually  and  simultaneously. 


Figures  2.4a,b,c  and  Table  2.2  give  typical  results.  Figure  2.4a  shows  a 
typical  sample  of  a  bi-spectral  ocean  background  histogram  -  a  very  high  peak  in  the 
low  intensity  region.  Although  we  designate  this  histogram  as  typical,  it 
should  be  realised  that  the  constant  "peak"  may  have  very  different  slopes  from 
location  to  location  in  the  same  image  due  to  the  inherent  inhomogeneity  of  the 
reflectivity  of  the  ocean  surface.  This  basic  problem  of  our  analysis  becomes 
even  clearer  when  one  looks  at  Figure  2,4b  where  we  show  the  histogram  of  the 
ocean  surface  slightly  obscured  by  a  light  haze.  The  features  of  the  surface 
histogram  are  still  there  but  one  already  sees  the  extension  with  a  definite 
slope  due  to  the  haze.  Finally,  Figure  2.4c  shows  the  full  histogram  of  a 


cloudy  region.  The  background  peak  Is  there  as  well  as  a  long  directed 
(correlated)  slope  culminating  in  a  very  high  peak  at  the  right  edge  of  the 
diagram,  both  due  to  the  presence  of  clouds.  It  is  clear  therefore  that  the 
problem  of  discrimination  between  clouds  and  background  is  a  major  one 
especially  for  the  smallest  clouds  with  varying  and  variable  optical  depths. 
Luckily,  the  analysis  is  not  too  sensitive  to  the  exact  value  of  the  threshold 
chosen.  One  last  method  was  also  tried  and  that  was  the  use  of  a  Laplaclan 
to  find  local  changes  in  the  brightness  gradient  between  a  given  pixel  and 
its  neighbours.  The  results  are  shown  in  Figures  2.5(a)  and  2.5(b)  for 
MSS4  and  MSS7  respectively.  The  Laplaclan  is  shown  as  a  function  of  the 
grey  level.  Going  back  from  high  grey  levels  we  find  positive  extrema  at 
grey  level  values  of  68  for  MSS4  and  at  2  for  MSS7,  quite  similar  to  those 
found  by  the  incremental  method. 

2.5  Summary 

In  summary  then,  none  of  the  methods  is  good  enough  by  itself  and  man- 


computer  interactive  choices  are  necessary  for  each  separate  image. 
Usually,  method  2.1  has  been  adopted. 
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The  statistical  properties  of  single  clouds 

The  next  stage  after  the  definition  of  the  cloud-background  demarcation 
intensity  in  a  given  spectral  image  is  the  treatment  of  the  clouds  as 
aggregates  of  "cloud"  pixels. 

All  pixels  with  intensities  above  or  equal  to  the  limit  are  termed  "cloud 
pixels",  all  others  are  "background". 

The  next  step  to  be  described  is  the  definition  of  a  cloud.  Afterwards, 
we  shall  deal  with  the  location  and  counting  of  clouds,  cloud  areas,  perimeters, 
holes  in  clouds,  inter-cloud  distances  as  well  as  various  derived  quantities  - 
like  the  fractal  dimension  of  a  cloud  field,  the  ratio  between  cloud  area  and 
its  surrounding  and.  the  various  probability  distributions  of  these  and  other 
quantities. 

-ii',  ,■  - .  ,f.  ;  ,  ■>  •'r  • 


The  first  step  in  this  long  line  of  activities  is  the  definition  of  a  , 

file  (A)  that,  contains  1  for  each  pixel  of  cloud  intensity  and  0  for  each  pixel 

'  •  •-  : i i  .  .  i  !  n..  :•  •  •  •  ■  .  • .  • >  .  .y  -t: 

of  background  .in teas ity..  }  in  .o^her  words  -  the  creation  of  a  high-contrast  map 
from  the  satellite  image.  ^ 

The,  neighbours,  of  a  given. pixel,  X..  in  row  i  and  column  j  in  a  cloud 

-  •  v!  ■  ■  i  .-xi  n  •,£  ,-r  -to !  r;«n  <  ■  •}.  r.iij.  <-«,;«•  n  >  .••;  ,  i  .. 

are  defined  as,  part  of  the  same-  cloi|d  if  t,hey  have  the  appropriate  intensity, 
viz,,  are  assigned  the  value  1.  ...  _ 

We  take  into  account  all  eight  neighbours  of  each  cloud  pixel  X...  ,  ,-1 

(Xi-l,j-r  Xi-1,j*  Xi-1  ,j+l  ’  Xi  ,j-l  ’  Xi,  j+1’  xi  +  i,  j-r  xi  +  i,j’  xi  +  i,  j  +  i} 
as  potential  parts  of  the  same  cloud. 
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3 . 1  The  enumeration_of  clouds 

Our  afm  here  (s  to  count  the  clouds  and  to  evaluate  their  area  In  pixels. 
At  the  end  of  the  process  to  be  described,  one  Identifying  cloud  number  will 
be  assigned  to  all  pixels  that  belong  to  the  same  cloud. 


These  numbers  -  the  "cloud  numbers"  -  replace  the  l's  on  file  A.  In 

addition  we  file  the  cloud  areas  as  a  vector. 

The  program  creates  a  fictitious  line  of  zeros  above  the  upper  edge  of  file  A. 

At  every  stage,  the  program  reads  one  line  Cm  }  of  l's  and  0's  from  file  A.  It 

*■  »  J 

always  uses  two  lines,  {M  }  and  {M  .}  simultaneously  in  the  counting 
pro cess. 


When  the  process  is  finished  for  one  line,  a  line  of  cloud  numbers  of 

{M. .}  replaces  {M  }  on  file  A.  Then  the  next  line  of  M  .  is  read  and 
ij  i»j  i,j+i 

the  process  repeated. 

Before  initiation  of  the  counting  process,  a  vector  JJ  is  zeroed.  |n  JJ., 
we  then  keep  the  cloud  number  connected  to  cloud  i.  Cloud  i  itself  may  be  t 

connected  to  another  cloud  number,  etc. 

For  every  line,  four  neighbours  are  defined  for  each  pixel  -  its  left 
neighbour  on  the  same  line  and  its  three  closest  neighbours  in  the  line  above. 

Pixels  on  the  left  margin  have  only  two  neighbours,  those  on  a  right  margin 
have  three  neighbours. 


.V; 


For  each  pixel ,  we  find  the  maximum  cloud  number  (MAX)  of  its  neighbours  and 


the  minimum  cloud  number  (KIN),  without  counting  the  0  (background)  pixels. 
Two  cases  may  occur: 


(a) 


Max  “0;  in  this  case  we  introduce 

on  line  {M..},  a  new  cloud  number 
U 

"cloud  number"  N  namely  ICL(N); 


into  the  location  of  the  pixel 
N,  and  enlarge  by  1  the  area  of 


(b)  Max  4  0;  we  have  to  check  whether  "cloud  number"  MAX  is  connected  to 
another  cloud  number.  If 


I  CL  (MAX)  -  0, 

MAX  is  infact  connected  to  another  cloud  number.  This  is  because  for  any  cloud 
number  that  is  found  to  be  connected  to  cloud  MAX,  ICL(MAX)  is  added  to  the  area 

of  that  cloud  and  is  then  put  to  zero.  Therefore  in  the  case 

I  CL (MAX)  =  0, 

MAX  is  exchanged  for  JJ(MAX)  and  we  check  if 

I  CL  (J J (MAX) )  =  0 

etc.  until  one  finds  a  cloud  number  for  which 

I  CL  (JJ(MAX))  4  0 

This  JJ  is  the  initial  cloud  number,  to  which  pixel  with  cloud  number  MAX  is 
connected.  MAX  is  replaced  by  this  JJ.  If  abinitio  ICL(MAX)  4  0,  it  is  not 
necessary  to  find  a  new  MAX. 


Now  again  two  different  cases  may  be  distinguished: 

(a)  The  new  MAX  cloud  number  equals  MIN  or  in  other  words  all  the  neighbours 
have  the  same  "cloud  number".  Then  we  introduce  the  new  MAX  into  the  location 

of  the  pixel  checked  out  on  line  {Mjj}  and  increase  the  area  of  the  cloud 

I CL (MAX)  by  one. 

(b)  The  new  MAX  is  different  from  MIN.  In  this  case,  one  has 

to  check  for  each  neighbour  if  its  cloud  number  is  connected  to  another 
cloud  number.  We  carry  out  the  above  process  until 

the  initial  cloud  number  is  found  to  which  the  neighbour  is  connected.  We  connect 
the  area  ICL  (initial  cloud  number  connected  to  neighbour)  to  the  area  ICL  (new 
MAX),  zero  ICL  (cloud  number  of  neighbour)  and  ICL  (initial  cloud  number  connected 
to  neighbour),  enter  the  new  MAX  into  JJ  (cloud  number  of  neighbour)  and  into 
JJ  (initial  cloud  number  connected  to  neighbour).  Then  one  enters  the  new  MAX 
into  the  location  of  the  pixel  checked  out  on  line  and  increases  the 

area  of  ICL(MAX)  by  one. 

When  the  process  is  finished  for  each  line,  one  writes  M. .  on  file  A  and 

i  J 

checks  for  value  of  maximum  cloud  number.  |f  the  largest  cloud  number  reaches 
within  500  of  the  dimension  of  the  vectors  JJ  and  ICL,  one  executes  a  reduction 
of  the  cloud  numbers  -  eliminating  "empty"  numbers  in  the  following  manner. 

We  first  find  the  Initial  cloud  number  to  which  every  cloud  number  is 
connected;  we  scan  JJ.  for  i  from  1  to  NMAX,  where  NMAX  is  the  largest  cloud 
number  reached  at  the  given  stage,  and  check  whether  JJ.  =  0.  If  not,  one  checks 
whether  JJ(JJ.)  is  equal  to  zero  and  so  on  until  equality  to  zero  is  reached. 
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The  index  of  JJ  for  which  the  equality  to  zero  is  reached  is  the  initial 
cloud  number.  One  writes  the  initial  ctoud  number  for  every  cloud  on  an 
auxi 1 iary  fi 1e. 

Then  one  scans  I  CL.  for  i  from  1  to  NMAX  and  finds  the  number  of  cases  of 

i 

ICL^  *  0  up  to  each  case  of  I  CL ^  t  0,  and  stores  it  in  ICLO^.  For  each  i  for  which 

ICL^  ^0,  we  write  on  an  additional  temporary  auxiliary  file  the 

value  of  ICLO^.  The  latter  quantity  is  the  number  that  has  to  be  subtracted 

in  order  to  eliminate  the  intervals  formed  by  the  cloud  numbers  for  which  ICL.  =  0. 

For  each  i  for  which  ICL.  ■  0,  we  write  on  this  auxiliary  file,  the  cloud 

number  i. 

Then  we  read  the  two  auxiliary  files  and  for  each  cloud  number,  i,  between 
1  and  NMAX,  we  subtract  from  its  cloud  number  the  appropriate  ICLO  « 

or  the  number  i  itself. 

The  result  is  a  vector  with  new  cloud  numbers,  arranged  sequentially  without 
gaps.  Part  of  the  new  numbers  are  zeros. 

We  read  the  file  A  once  more  upto  the  line  reached  in  the  counting  process 

and  enter  instead  of  each  M. .  the  new  cloud  numbers  and  rewrite  file  A. 

'  J 

Finally,  we  sort  ICL.  according  to  the  new  (reduced)  cloud  numbers.  We 
scan  ICL.  and  sort  from  i  =1  all  the  areas  ICL.  that  are  different  from  zero. 

i  i 

The  rest  of  the  vector  ICL.  is  set  equal  to  zero. 

i 
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At  the  end  of  the  reduction  process,  we  define  the  new  maximal 
cloud  number  as  NMAX,  check  if  the  last  line  of  file  A  is  reached, 
read  the  last  processed  line  of  file  A  to  M  ,,  read  another  line 
from  file  A  and  continue. 

A  flow  chart  of  the  complicated  process  described  above  is  given 
in  Figure  2.6. 
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3.2  Holes  in  clouds 

The  identification,  counting  and  filing  of  holes  is  analogous 
to  the  operation  with  cloud  areas.  In  order  to  forestall  the  possibility 
of  pixels  in  a  given  cloud  that  are  empty  but  cannot  be  defined  as  part 
of  a  hole,  we  count  here  as  neighbours  only  those  pixels  that  have  a  common 
boundary  (namely  for  pixel  X.^.  -  X.+1>j,  Xj-i  jXj,j-i'  Xi,j+I^' 

Using  the  file  of  cloud  numbers,  we  identify 
those  pixels  inside  a  cloud  that  are  not  cloud  pixels. 

We  end  up,  with  the  number  of  holes  and  their  association  to  clouds. 
Figure  2.7  shows  typical  results  (Case  7032).  The  size  distribution  of  the 
holes  has,  similarly  to  that  of  the  clouds,  a  power-law  type  of  behaviour. 

The  slope  is  however,  much  steeper.  In  the  case  shown,  the  slope  of  the 
size  distribution  of  the  holes  is  12,  versus  a  value  of  2.7  for  that  of  the 
cloud  sizes.  The  influence  of  the  holes  on  the  results  of  our  analysis  should 
therefore  be  minor,  except  for  the  very  smallest  clouds.  It  is  however  always 


taken  into  account. 
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3.3  Cloud  perimeters 

Cloud  perimeters  are  determined  for  each  cloud  by  adding  the  number  of 
outside  edges  of  pixels  in  the  x  direction  times  their  length  Ax,  to  the 


number  of  such  edges  in  the  y  direction  times  the  pixel  length  in  the 
y-direction,  Ay.  So  the  i  th  cloud,  has  a  perimeter  P.  given  by: 


n  m 

P.  =  (I  P)  Ax  +  /IP) Ay 
1  ix=l  X  iy=l  ^ 


All  edges  abutting  non-cloud  areas , incl uding  holes,  are  counted. 

We  do  not  take  into  account  that  non-cloud  pixels  abutting  on  cloud  pixels 
may  be  partially  cloudy.  This  determination  of  possible  partial  cloudiness  of 
such  single  pixels,  which  are  below  our  resolution,  is  impossible,  except  by 
some  assumption  -  e.g.  50%  coverage.  Such  an  assumption  does  not  increase  our 
accuracy  and  may  lead  to  an  under-estimate  of  both  perimeter  and  area  on  our 
part.  On  the  other  hand,  the  fact  that  we  do  not  smooth  the  serrated 
cloud  edges  leads  to  a  possible  over-estimate  of  perimeter  and  area  of  about 
the  same  magnitude. 

As  the  discrimination  between  cloud  and  non-cloud  pixels  is  indeterminate 
to  about  one  or  two  gray  levels,  we  do  not  think  it  fruitful  to  make  any 
correction  for  the  above  effects. 
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3.4  Cloud  areas 

The  area,  A,  of  a  cloud  I,  Is  simply  defined  as  the  number  of  pixels 
connected  to  it  multiplied  by  the  area  of  a  single  pixel,  viz: 


A.  -  (  E  a.)  Ax  Ay 
1  i-1  ' 


3.5  Cloud  centers 

The  center  of  a  cloud  is  defined  as  its  "center  of  gravity",  where  equal 
weights  are  given  to  all  pixels.  Therefore,  the  x  coordinate  of  the  cloud 
cen ter  is: 


C  «  E  x,/N 


and  similarly,  the  y-coordinate  is 


C  =  E  y./N 
y  i-1  ' 


*  y s v  !->’av,!vw.v  v  *  -v 
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3.6  Minimal  distance  between  cloud  pairs 


The  calculation  of  minimal  inter  cloud  distance  is  based  on  the  assumption 


that  each  cloud  is  defined  by  its  center  coordinates.  Therefore  the  distance 


between  two  clouds  is 


0  ‘  <  <Cx2  -  Cx,)2  +  (Cy2  -  Cyl>2) 


Each  cloud  has  at  least  one  neighbour  that  is  closer  to  it  than  all  the  others. 


The  distance  between  the  cloud  in  question  and  its  closest  neighbour  defines  the 


linimal  inter-cloud  distance,  DISMIN,  for  that  cloud. 


Then,  we  define  a  second  minimal  distance  D  ,  namely  the  minimal  distance 


between  a  cloud,  i,  and  its  neighbours,  j,  of  _>  1.1  . 


This  is  to  take  into  account  that  a  given  cloud  may  have  very  small  close 


neighbours  in  the  image  that  may  physically  be  part  of  the  larger  cloud. 

3.7  The  ratio  of  cloud  area  to  cloud  area  of  influence 

For  the  latter  distribution  of  minimal  distances,  we  calculate  for  each 


cloud  a  relationship  between  the  cloud  area  nr  ,  and  an  associated  cloud  - 

e 


f  ree  "r  i  ng"  around  it 


(Dmin)^  2 
-= —  -  r 
2  e 


Vfe  assume  that  one  half  of  the  minimal  distance  between  the  cloud  and  its 


neighbour  defines  the  cloud  area  of  i nf 1 uence, and  that  clouds  are  circular  of 


area  A, 


A  =  irr 


Oc.j 


These  ratios  must  of  course  be  positive.  Let  us  investigate  the 


possibilities  for  negative  values* 


The  ratio  y  is  equal  to  zero  only  f°r  re  “  ®  (no  cloud) 


Negative  values  y  <  0  occur  when 

(Dmin)2  ~r  2  <  0 
2  e 


or 


Dmin  „ 

2  e 


This  may  occur  when  two  elongated  clouds  occur  side  by  side. 

In  all  cases  analysed,  almost  all  data  show  positive  values  which  are 
almost  independent  of  cloud  size. 
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he  composite  histogram  method.  The 
evel  threshold  value.  The  Hi,  H2  an 
.n  the  text. 


Figure  2.3: 


The  effect  of  noise  on  Hi,  H2  and  H3  and  on  the  chosen  threshold 
rfTP,  HIM,  H3P  and  H3M  are  elaborated  upon  in  the  text. 
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Figure  2.  A: 


FLOW  CHART  OF  THE  ENUMERATION  OF  CLOUDS 


Figure  2.7:  The  size  distribution  c  holes  in  clouds. 
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Table  2.1:  Variation  of. cloud  parameters  with  grey  level  threshold 
Image  no  .  7025,  ljnes  600  -  900,  columns  1300  «-  1800 


Threshold 
grey  level 

No.  of 
clouds 

No.  of 
small  clouds 

No.  of 
large  clouds 

Cl oudiness 
in  % 

Cloudiness 
larqe  clouds 

MSS  -  4: 

Threshold  Chosen 

-  69 

80 

511 

146 

365 

13.2 

13.0 

78 

512 

143 

369 

13.7 

13.5 

76 

513 

144 

369 

14.2 

14.0 

74 

514 

134 

380 

14.7 

14.5 

72 

522 

126 

396 

15.3 

15.1 

70 

511 

122 

389 

15.9 

15.8 

68 

514 

124 

390 

16.6 

16.4 

66 

511 

122 

389 

17.3 

17.2 

64 

513 

124 

389 

18.2 

18.0 

62 

528 

129 

399 

19.1 

19.0 

60 

559 

165 

394 

20.3 

20.1 

MSS  -  7: 

Threshold  Chosen 

-  2 

10 

405 

130 

275 

5.8 

5.8 

9 

411 

121 

290 

6.4 

6.2 

8 

419 

112 

307 

7.1 

7.0 

7 

443 

118 

325 

7.9 

7.7 

6 

449 

115 

334 

8.8 

8.7 

5 

479 

135 

344 

10.0 

9.9 

4 

500 

148 

352 

11.4 

11.2 

3 

514 

135 

379 

13.2 

13.0 

2 

544 

149 

395 

15.8 

15.6 

1 

1961 

546 

23.5 

21.6 

: 
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TABLE  2.2:  Numerical  values  of  bi-spectral  histograms  of  Figure  4 
A.  PART  OF  IMAGE  7025  CONTAINING  CLEAR  OCEAN 
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Chapter  3»  The  distribution  of  Cumulus  cloud  sizes  and  intercloud 
distances  as  derived  from  LANDSAT _ 

3-1  Introduction 

It  is  widely  assumed  that  clouds  act  as  important  moderators  of  the 
Earth's  radiation  balance  on  both  regional  and  global  scales.  Solar  and 
thermal  radiation  fields  in  the  atmosphere  are  significantly  modified  by 
the  presence  of  cloud  fields  (e.g.  Welch  et  al,  1980,  Davies  1978,  Joseph 
and  Davies  1985,  Harshvardhan  and  Thomrs  1984). 

The  radiational  effects  of  individual  clouds  depend  on  many  factors  - 
cloud  thickness,  cloud  top  height,  cloud  vertical  and  horizontal  geometry, 
the  illumination  geometry  as  well  as  liquid  water  content,  cloud  droplet 
size  distribution  and  atmospheric  water  vapor  profiles.  For  cloud  amounts 
larger  than  about  30  percent,  several  types  of  cloud-cloud  interactions 
due  to  the  "brokenness"  of  the  cloudiness  become  important  (Welch  et  al, 

1980,  Joseph  and  Davies  1985,  Gube  et  al  1980,  Harshvardhan  and  Weinman 
1982,  Harshvardhan  and  Thomas  1984). 

Clouds  may  also  be  used  as  diagnostic  features  in  studies  of  weather 
and  climate.  The  cloud  fields  affect  temperature,  humidity  and  motion 
patterns  of  their  ambient  atmosphere  (e.g.  Paltridge  1974a, b,  Cox  et  al,  1973, 
Reynolds  and  Cox  1975,  Stephens  1976,  Fravalo  et  al,  1981,  Schubert  et  al, 
1977,  Yanai  et  al,  1977,  Foltz  and  Gray  1979,  Cho  1982,  Gray  and  Jacobson 
1977,  Fingerhut  1977). 

Therefore  any  study  of  atmospheric  circulations  and  their  accompanying 
phenomena  must  take  into  account  the  cloud-fields  both  as  to  their  role  in  the 
dynamics  of  the  development  and  their  role  as  diagnostic  features. 


Major  efforts  on  the  definition  and  parameterization  of  cloudiness  are 
underway  (e.g.  ICCPP) .  We  are  engaged  in  a  detailed  study  of  the  morphology 
of  broken  cloud  fields  as  viewed  with  high  resolution  from  satellites.  The 
morphology  of  a  cloud  field  is  the  geometrical  structure  of  the  clouds  and 
their  spatial  inter-relationships.  As  stated  before,  the  morphology  of 
cloud  fields  is  strongly  connected  to  atmospheric  thermo-hydrodynamics,  there¬ 
fore  cloud  and  cloud  field  properties  can  serve  as  both  diagnostics  and  pre¬ 
dictors  of  atmospheric  phenomena.  However,  the  morphology  of  clouds  and 
cloud  fields  also  has  important  effects  on  radiative  transfer  through  a 
planetary  atmosphere.  These  are: 

(a)  Radiative  transfer  of  energy  fluxes  -  e.g.  heating  phenomena, 
availability  of  solar  energy  at  the  Earth's  surface; 

(b)  Radiative  transfer  of  information  -  e.g.  detection  of  surface 
features  and  properties  through  cloudy  layers. 

In  this  paper  we  summarize  some  of  the  initial  results  of  our  ongoing 
analysis  of  high  resolution  (70x80  m)  of  LANDSAT  data  taken  over  fair  weather 
cumulus  cloud  fields  in  different  parts  of  the  world. 
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LANDSAT  and  similar  high  resolution  data  have  a  special  place  in  the  study  of 
cloud  field  features.  Detailed  analysis  of  the  high  spatial  frequency  variations 
of  the  cloud  field  in  parallel  with  consideration  of  low  frequency  (e.g.  meteoro¬ 
logical  satellite  data  will  make  it  possible  to  take  into  account  the  high  frequency 
effects  of  cloud  fields  knowing  only  the  meteorological  satellite  data.  In  addition 
studies  of  convection  and  planetary  boundary  layer  dynamics  will  be  much  supported 
by  high  resolution  satellite  data. 

The  imagery  available  to  us  is  in  the  form  of  CCT's  with  a  data  density  of 
6250  b.p.i.  of  all  spectral  bands  as  well  as  of  black  and  white  pictures.  Imagery 
of  suitable  quality,  highest  in  all  bands,  and  cloudiness  0.1  to  0.6  is  chosen  from 
EROS  catalogues,  principally  over  the  oceans  at  this  stage  of  our  study. 

After  reading,  copying  and  checking  the  quality  of  the  tapes,  we  Identify  from 
the  black  and  white  pictures  of  two  bands  (4  and  7  in  the  case  of  LANDSAT  S2,  3), 
the  location  of  the  cloudy  and  clear  areas  to  be  analysed. 

Histograms  of  the  intensity  from  these  areas  are  then  constructed  from  which 
we  determine  first  whether  the  intensities  of  the  cloud  images  are  saturated  or  not, 
viz  whether  there  is  a  large  peak  at  the  highest  gray  levels.  If,  as  is  often  the 
case,  the  cloud  images  are  in  fact  saturated,  we  determine  the  lower  limit  of  cloud 
pixel  intensities  in  the  following  way.  The  lower  limit  of  cloud  pixel  intensity 
is  iteratively  reduced  from  a  value  close  to  the  maximum  in  both  bands  till  the 
following  conditions  are  met: 

(a)  the  total  cloudiness  approaches  the  same  limit  in  the  two  bands, 

(b)  the  total  number  of  clouds  in  both  bands  is  approximately  equal, 

(c)  the  ratio  of  the  number  of  "small  clouds"  (less  than  four  pixels  in  area) 

to  total  number  of  clouds  increases  suddenly,  viz  the  noise  level  has  been 
reached. 

The  lower  limit  varies  from  image  to  image  and  must  be  arrived  at  in  the 
above  manner. 


If  the  clouds” are  not  saturated,  we  use  one-  and  two-dimensional 
histograms  of  cloudy  and  cloud-free  areas  to  define  the  initial  cut  between 
clouds  and  background  and  then  do  the  analysis  described  above.  Additionally, 
if  the  intensities  are  not  saturated  we  use  data  from  bands  4  and  6  (in  the 
case  of  Landsat  3)  or  equivalent  bands  in  Landsats  4  and  5  to  determine  the 
cloud  height  in  the  following  manner. 

The  height  of  cloud  tops  is  determined  by  assuming  that  the  cloud 
reflectivity  does  not  vary  spectrally.  The  cloud  reflectivity  can  then  be 
deduced  from  the  red  or  solar  I.R.  band  and  used  in  the  blue  band  to  calculate 
the  Rayleigh  optical  depth  of  the  atmosphere  above  the  cloud.  From  the 
optical  depth,  the  pressure  and  the  height  are  then  deduced.  A  separate  paper 
will  be  published  comparing  this  method  with  others  -  e.g.  use  of  thermal  I.R- 
intensities  in  LANDSATS  4  and  5  (Joseph,  Mekler,  Kahn  and  Podolak,  1986). 

The  noise  of  the  image  was  analysed  and  we  made  several  tests  on  the 
efficacy  of  smoothing  e.g.  de-spiking,  but  the  effect  turned  out  not  to  be 
important  for  our  purposes  and  consequently  was  deleted  from  our  procedures. 

The  process  is  summarized  in  Figure  3.1. 
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The  f irstr  type  of  clouds  to  he  studied  is  Fair  Weather  Cumulus  over 
sea  in  low  and  sub-tropical  latitudes. 

The  reasons  for  this  choice  are  several.  Firsts  it  is  easiest  to  develop 
an  analysis  technique  for  small,  discrete  and  bright  objects  over  a  darker 
and  usually  fairly  reasonably  homogeneous  background.  This  type  of  clouds  is 
therefore  of  intrinsic  interest  as  well  as  a  separate  species  of  clouds. 

Our  plan  is  to  extend  the  analysis  to  clouds  over  land  as  well  as  to 
include  pattern  recognition  of  extended  meso-  and  macro-scale  cloud  regimes. 

In  addition  to  the  basic  single  cloud  properties  described  in  this  study  - 
size  and  inter-cloud  distance  distributions  and  total  cloudiness  -  we  derive  also 
cloud  perimeters  and  cloud  top  heights  as  described  above,  and  the  distribu¬ 
tion  of  holes  in  clouds  (see  next  chapter) . 

In  addition,  we  derive  multi-variate  properties  -  the  fractal  dimension 
of  the  cloud  field,  the  relations  between  the  size  and  inter-cloud  distance, 
the  contribution  of  clouds  of  different  sizes  to  total  cloudiness  for  16 
intensity  ranges,  the  ratio  of  cloud  to  cloud  surrounding  areas  as  a  function  of 
cloud  areas,  etc.  (see  next  chapter). 

Finally,  we  connect  the  various  statistics  to  atmospheric  dynamics  in 
several  ways. 

All  these  will  be  presented  in  subsequent  papers  for  a  variety  of  locations 


covering  the  globe. 
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3.3  The  Size  Distribution  of  Clouds 

Typical  data  on  the  experimental  size  distribution  are  given  in  Figure  3,3 
Initially,  we  fitted- the  data  to  an  exponential  curve  in  order  to  compare  our 
results  to  previous  work  (e.g.  Plank  1969).  The  part  that  was  fitted  is  from 
the  minimum  size  -  four  pixels  or  r  min  ^  64  m  -  to  a  maximum  size  prescribed 
by  that  effective  maximum  radius  beyond  which  the  data  degenerate  to  an  irregularly 
distributed  tail  of  single  clouds. 

It  is  clear  from  the  figure  that  a  power  law  describes  the  data  better  than 
an  exponential  one.  In  addition  it  should  be  mentioned  here  that  a  power  law 
is  concomitant  with  the  fractal  nature  of  clouds  (Mandelbrot,  1982)  which  will 
be  dealt  with  in  a  subsequent  paper.  (Cahalan  and  Joseph  1986). 

All  our  results  on  the  cloud  radius  distributions  are  summarised  in 
Table  3.1  and  in  the  Appendices  to  this  report.  This  table  includes  the  data 
from  previous  studies  of  U-2  photographs  (Plank  1969)  as  well  as  from  that  of 
Hozumi  (1982)  for  cloudiness  less  than  20%.  The  other  data  from  Hozumi  are 
for  larger  amounts  of  cloudiness  in  which  the  type  of  clouds  might  be  different 
from  ours  -  Fair  Weather  Cumulus. 

In  order  to  get  a  visual  impression  of  the  data  and  the  possible  range  of 
variability  we  present  comparisons  of  the  size-distribution  and  the  accumulated 
size  distribution  (integrated  backwards)  for  three  cases  spanning  the  range  of 
cloudiness  in  our  study  -  0.5  to  about  40  percent  in  figures  3.3(a)  to  3.3(j). 

For  each  band  two  figures  are  shown  -  the  size  distribution  and  the  accumulated 
size  distribution  (integrated  backwards  from  the  largest  cloud).  The  actual 
data  points  are  situated  at  the  top  of  each  L-shaped  segment  of  the  size- 
distribution  histogram  and  at  the  right  hand  edge  of  the  base  of  the  "L"  shape 
in  the  case  of  the  accumulated  distribution  curve.  The  fit  to  the  size  distri¬ 
bution  is  indicated  by  a  sloping  line  and  the  cut-off  point  by  a  vertical  line 
intersecting  the  radius  coordinate.  The  largest  clouds  may  not  always  appear 
in  the  size  distribution  but  do  appear  in  the  accumulated  size  distribution. 

For  the  sake  of  comparison  all  figures  are  with  identical  scale.  In  general 
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it  is  clear  from  the  three  figures  as  well  as  from  the  Appendices  that  the 
size  distributions  are  quite  similar  for  all  amounts  of  cloudiness  larger 
than  about  two  percent  .  In  addition,  it  is  evident  that  for  each  image,  the 
two  bands  give  the  same  result.  This  constitutes  one  of  the  checks  in  the 
analysis. 

The  slopes,  beta,  of  the  three  size  distributions  shown,  with  respect 
to  the  effective  radius,  are  almost  equal  -  the  average  value  of  the  slope 
of  all  our  data  is  2.83  ±  0.20  taking  into  account  ours  as  well  as  Plank's 
relevant  data  (Hozumi's  data  were  not  taken  into  account  because  they  describe 
mixed  clouds).  This  is  obvious  from  looking  at  the  Figures  3.3(a)  to  (j)  and 
the  Tables.  Only  the  integration  constants,  K,  differ.  The  value  of  the 
slope  is  larger  than  the  value  of  7/3  which  follows  from  the  clouds  being  of 
fractal  nature  with  a  dimension  of  the  perimeter  of  4/3  (Lovejoy  1982). 

Maximal  cloud  size  varies  between  0.3  and  2.0  km. 

The  constant  K  varies  much  more  from  location  to  location.  The  present 
data  show  no  strong  correlation  of  K  with  P  -  the  percentage  of  cloudiness  - 
or  with  any  other  parameter. 

It  must  be  especially  emphasized  that  the  slope  w.r.t.  the  areas.  A,  in 
a  log-log  diagram  is  half  the  previous  value  -  namely  about  1.43.  This  is 
important  and  relevant  when  one  tries  to  fit  the  present  results  to  a  cloud 
formation  model  by  convection  or  by  aggregation  based  on  a  fractal  approach  - 
e.g.  D.L.A.  of  clusters  (Family,  Private  Communications , 1986) .  It  seems  that 
the  above  results  as  well  as  the  fractal  dimensions  of  the  perimeter-in  our  case 
about  1.3  and  1.55  as  shown  by  us  in  the  next  chapter-make  it  reasonable  to 
expect  a  good  fit  to  one  of  the  cluster  aggregation  models,  (ibid,  Private 
Communication) . 
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Another  interesting  feature  is  that  the  maximal  cloud  radius  in  each 


case  seems  to  be -correlated  to  the  cloud  amount  -  the  correlation  coefficient 


is  over  0.95.  For  the  higher  cloud  amounts,  it  is  possible  that  the  largest 


clouds  (larger  than  10  km  )  may  have  a  different  frequency  distribution  of 


the  cloud  areas.  This  is  apparent  only  in  the  accumulated  size  distributions 


given  in  the  figures  -  viz.  there  may  be  a  change  in  slope  beyond  10  km 


3.4  The  Distribution  of  Minimal  Distances  Between  Pairs  of  Cloud  Centers 


For  each  cloud  in  a  given  sample,  we  find  the  coordinates  of  its 


"center  of  gravity".  Then,  for  any  given  cloud,  we  find  the  one  other  cloud 


that  is  closest  to  it  in  terms  of  the  distance  between  the  two  "centers  of 


gravity"  in  two  circumstances  -  first  all  clouds  are  taken  into  account  and 


secondly  only  those  clouds  that  are  at  least  one  tenth  larger  than  the  area  of  the 


cloud  in  question.  One  example  of  the  distribution  of  these  distances  is  given 
in  Fig.  3.4a  and  3.4b  in  the  two  spectral  bands  analysed,  MSS 4  and  MSS7 , 


showing  again  that  the  same  results  are  obtained  in  both. 


The  distribution  is  strongly  peaked  -  in  both  spectral  bands  analysed  -  at 


a  value  of  about  0.5  km.  This  peakedness  is  a  feature  common  to  all  distributions. 


Some  representative  data  are  given  in  Table  3  in  the  form  of  fits  to  log-normal 


distributions.  This  seemed  to  be  a  natural  choice  due  to  both  the  shape  of  the 


curves  as  well  as  to  the  common  occurrence  of  that  distribution  in  many  studies 


of  cloud  physics  (e.g.  Lopez  1977). 


In  addition,  the  log-normal  distribution  may  be  a  result  of  the  cloud 


formation  process  (e.g.  Lopez  1977) 


A  fairly  surprising  result  is  that  the  values  of  PI  and  (especially)  PI 


are  very  similar  for  all  locations  -  not  only  those  shown  in  the  table.  Thi 


may  be  due  to  similar  conditions  of  convection  at  all  locations  -  mid  mornir 


fair  weather  cumulus. 
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The  fit  to  a  log-normal  may  however  not  be  the  optimal  approach.  In  a 
separate  note,  we  show  (Cahalan  and  Joseph  1986)  that  for  a  random  distribution 
of  cloud  centers  in  the  plane  one  should  expect  a  distribution  of  a  different 
type.  We  therefore  decided  to  try  and  fit  the  experimental  histogram  of  the 
distances  between  nearest  neighbours  using  different  assumptions.  We  assumed 
that  the  probability  of  finding  a  cloud  center  at  distance  r  from  any  other 
one  is  given  by 


W(s)  ds  *  [1  -  W(s')  ds’]  2TTsn  f(s)  ds 


where  the  expression  in  square  brackets  gives  the  probability  of  not  finding 


a  cloud  at  a  distance 


s'  <  s 


27T  snds 


gives  the  random  probability  of  finding  a  cloud  center  in  an  annulus  of  radius 
s  of  area  2tt  s  ds  and  when  the  average  density  of  cloud  center  is  n.  The 
function  f(s)  describes  the  deviation  from  a  random  distribution  in  the  plane. 

Equation  (l)  is  solved  by  us  and  its  implications  are  more  fully  discussed 


in  another  paper  (Cahalan  and  Joseph  1986).  Here  we  present  results  assuming 


f(s)  =  s 


(N—  1 ) 


V/f  »  S' 


L.JS  . 


and  that  therefore 


W(s)  *  2nn  sN  exp{-2ir  ns^*"^ /(N+l)} 

Therefore  N  >  1  implies  deviation  from  a  random  distribution  of  clouds  in 
the  plane,  such  that  the  probability  of  finding  a  cloud  distanced  s  from 
the  center  cloud  increases  with  s.  We  therefore  show  the  same  data  as 
before  but  fitted  to  such  a  function  in  Table  4. 

Again  the  similarity  of  n  and  b  at  different  locations  are  often 

striking,  but  in  addition  the  fact  that  n  is  often  close  to  one  shows  that 

the  distribution  of  cloud  centers  in  the  plane  is  then  close  to  random. 

The  reason  for  the  latter  conclusion  is  that  a  purely  random  distribution 

would  result  in  a  value  of  one  for  the  power  N  (Cahalan  and  Joseph  1966). 
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In  order  to -give  a  further  impression  of  the  features  of  the  distribution 
of  minimal  distances  between  cloud  centers,  we  present  Figures  3.5(a),  (b)  and 
(c)  for  amounts  of  cloudiness  of  1.4,  13.5  and  32.9  percent  respectively.  The 
figures  are  for  one  band  only  (MSS4)  as  the  results  for  the  other  one  are 
similar.  The  histograms  are  not  normalised  but  it  is  easy  to  check  that 
normalisation  to  a  constant  number  of  "large  clouds"  (cloud  with  areas  larger 
than  four  pixels),  will  make  all  three  peaks  coincide  within  the  experimental 
error. 

Superposed  on  the  experimental  data  -  the  crosses  -  are  the  two  fits  that 
we  used  -  the  log  -  normal  and  the  Weibull  distributions.  The  two  fits  give 
similar  results,  except  for  the  region  of  minimal  distances  larger  than  2  km, 
where  the  Weibull  function  decreases  too  rapidly.  One  may  conclude  that  the 
log-normal  distribution  is  the  one  to  adopt,  similarly  to  the  usage  in  most 
other  studies  on  cloud  properties  (Lopez  1977).  However,  in  our  opinion, 
the  experimental  data  are  not  good  enough  for  such  a  conclusion  -  they  indicate 
»  a  distribution  with  a  given  peak  and  width  only.  It  may  however,  be  important 

to  reach  a  definite  conclusion  because  of  the  physical  phenomena  behind  each 
distribution  (e.g.  ibid  1977). 

We  investigated  statistically  the  possible  correlations  between  the 
amount  of  cloudiness,  the  number  of  "large"  and  "small"  clouds  and  their  ratio 
and  the  coefficients  of  the  two  fits  (See  Appendices). 


The  only  significant  correlation  is  that  between  the  amount  of  cloudiness 
and  N,  the  exponent  of  the  Weibull  distribution.  It  is  low  when  one  takes 
into  account  all  data  -  0.31.  When  one  successively  increases  the  lowest 
amount  of  cloudiness  taken  into  account,  the  correlation  increases  rapidly 
and  reaches  0.87  for  cloudiness  larger  than  15  percent  and  0.92  for  cloudiness 
larger  than  20  percent.  The  values  of  cloudiness  for  which  N  is  equal  to  one 
an  indication  of  randomness  in  the  distribution  of  cloud  centers  -  varies 
between  10  and  18  percent.  This  is  reasonable  in  terms  of  experience.  Fair 
Weather  Cumulus  clouds  usually  exhibit  a  random  distribution  for  such  low 
amounts  of  cloudiness. 


3.5  Discussion  and  Summary 

This  study  presents  the  initial  results  of  an  analysis  of  the  statistical 
distribution  functions  of  cloud  size  and  minimal  inter-cloud  distances  at  seventeen 
locations  spread  over  the  world's  oceans.  These  single  cloud  properties 
are  a  necessary  first  step  in  a  study  of  the  complex  morphology  of  broken  cloud 
fields  and  the  application  of  the  latter  to  relevant  problems  in  atmospheric 
physics. 

The  type  of  clouds  analysed  -  Fair  Weather  Cumulus  -  is  an  important  widely 
and  frequently  occurring  type  and  represents  the  prevalence  of  convection  in  the 
atmospheric  boundary  layer.  Fair  Weather  Cumulus  clouds  usually  have  aspect  ratios - 
of  about  one  (e.g.  Flank  1969,  Gagin  1983  private  communication)  so  that  our  cloud 
radii  distributions  are  also  representative  of  the  thicknesses  of  these  non- 
raining  clouds.  The  cloud  thickness  will  be  proportional  to  twice  the  effective 
radius  found. 

Some  interesting  commonalities  have  been  observed  in  the  distributions  found 
at  different  locations.  First,  a  power  law  fits  the  size  distributions  best. 

The  power,  B,  is  fairly  similar  from  location  to  location,  2.83  ±0.20,  and  is 
slightly  larger  than  7/3,  the  value  resulting  for  clouds  assumed  to  be  self-similar 
fractals . 

— s 

The  constant  K  of  the  power  law,  Kr  varies  more  from  location  to  location 
and  is  weakly  correlated  to  the  cloud  amount  as  shown  in  Section  3. 


The  number  of  clouds  tends  to  decrease  with  the  amount  of  cloudiness 
similarly  to  the  previous  results  from  different  times  of  day  at  the  same 
location  (Plank.  1969) .  This  brings  us  to  a  basic  deficiency  in  the  LANDSAT 
data,  namely  the  fact  that  they  are  available  at  one  time  of  day  only  -  mid- 
morning.  This  is  of  course  partly  off-set,  but  not  corrected,  by  the  possibility 
of  analyzing  many  locations  globally  using  the  same  instrumentation. 

In  the  future,  our  data  in  order  to  be  more  fully  useful,  will  need  to  be 
augmented  by  observations  at  other  times  of  day  (e.g.  by  use  of  airplanes). 

In  addition,  it  is  important  to  extend  the  lower  limit  of  resolution  by 
analysis  of  imagery  from  other  sources  (Thematic  Mapper,  SPOT  Satellite  and 
airplanes) .  This  point  will  be  gone  into  more  fully  in  subsequent  papers  dealing 
with  the  fractal  property  of  the  cloud  field  and  its  applications. 

The  fit  of  the  minimal  inter-cloud  distances  to  a  function  given  by: 

W(s)  =  2tt  nsN  exp,  <  -  Zn  ns^N+^  /  (N+l),  >  shows  that  in  many  cases 
studied  the  areal  distribution  of  clouds  is  close  to  random  -  N  ^  1, 
especially  for  cloudiness  lower  than  20  percent.  This  same  power,  N,  may  thus  serve 
as  an  indication  of  the  deviation  from  randomness  e.g.  the  cloud  field  west  of 
Israel  over  the  Mediterranean.  Much  additional  information  on  the  nature  of  the 
convection  may  be  derived  and  will  be  reported  on  separately. 

Further  analysis  of  imagery  at  other  locations  is  in  progress  as  well 
as  various  applications  to  problems  in  atmospheric  physics. 
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One  other  point  which  should  be  emphasized  here  is  the  extreme  difficulty 
of  developing  a  reliable  semi-automatic  data  analysis  scheme.  Great  effort 
went  into  the  latter's  development  and  a  significant  part  of  the  project  time 
was  spent  studying,  checking  and  implementing  the  various  possible  methods. 

Our  U.S.  counterpart  took  upon  himself,  inter  alia,  the  analysis  of 
different  aspects  of  the  problem,  e.g.  the  perimeter  -  area  relationships  of 
the  smallest  clouds,  close  to  the  limit  of  resolution  -  and  consequently  was 
able  to  analyse  in  depth  one  image.  His  study  was  however,  fundamental 

to  the  development  of  our  approach  in  Israel,  and  will  be  reported  separately. 
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FIGURE  3.1:  ANALYSIS  OF  LANDSAT  IMAGERY 
BASIC  PROBLEM  OF  DATA  ANALYSIS; 
DISCRIMINATION  BETWEEN  CLOUDS  AND  BACKGROUNDS 


BASIC  TOOLS 


Histogram"  of  Spectral  Intensity 
Bands  4,  6  and  7  used 
(generally  4  and  7  only) 


CLOUD 

INTENSITY 

SATURATED 


Scheme  I 

Lower  limit  of  cloud  pixel 
intensity  iteratively  re¬ 
duced  from  value  close  to 
maximum  intensity  in  each 
band  until: 

(a)  total  cloudiness  approaches 
identical  limit  in  the  two 
bands; 

(b)  ratio  of  "small-clouds" 
area  (less  than  4  pixels 

in  area)  to  total  cloud  area 
increases  suddenly  -  noise 
level  reaclTlL#. 


Scheme  II 

ID  and  2D  histograms  and 
spatial  coherence  of  Imagery 
used  together  with  similar 
histograms  of  selected 
background  areas  (e.g.  surface) 
to  define  cut  between  clouds 
and  backgrounds. 

Iterations  carried  out  till 
conditions  (a)  and  (b)  of 
Scheme  I  are  met 


Note:  Height  cannot  be  determined 
when  bands  saturated 


Size  distribution  and  accumulated  distributions  for  three 
Spanning  the  range  of  observed  cloudiness 
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Figure  3.4; 


Histograms  of  minimal  distances  between  cloud  centers  for 
hands  MSS4  and  7  in  case  7032  (See  Appendix). 
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The  distribution  function  of  minimal  distance  between  cloud 
centers  for  three  cases  spanning  the  range  of  cloudiness 
observed  -  (a)  Case  7067,  cloudiness  1.4%;  (b)  Case  7025, 
cloudiness  13.5%;  (c)  Case  7024,  cloudiness  32.9%. 
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TABLE  3.2:  The  size  distribution  of  Cumulus  clouds  over  oceans 


dN(re)  . 

d  r 

e 

K  r 

e 

where 

A  Is 

the  area  of  cloud 

Location  and  Center  Coordi. 

Band 

K 

8 

r  (km)  No. 

max 

of  clouds 

P(*)* 

A (km2)  ** 

Ocean  E.  of  Florida 

MSS4 

0.12 

2.86 

0.46 

2521 

32.9 

8.036 

N23°04,W078°44‘ 

MSS7 

0.11 

2.87 

0.46 

2420 

33.2 

Atlantic  Ocean  E.  of  Florida 

mm 

2.97 

0.46 

7329 

13.5 

17,172 

N24°30‘W79ol*7' 

mm 

EBB 

2.96 

0.46 

71*36 

12.9 

Atlantic  Ocean  E.  of  Florida 

331*9 

m 

17,235 

N24°30'W79°1*7 1 

mm 

3**  58 

BBS! 

Atlantic  Ocean  £.  of  Florida 

mm 

3.03 

0.53 

4738 

19.0 

15,200 

N24°30 ' U8 1 0 1 1 ' 

eH 

3.02 

0.43 

4671 

18.5 

Atlantic  Ocean  E.  of  Florida 

in 

mm 

3.01 

0.28 

1*235 

8.8 

16,243 

N2ito28'U79olt0' 

mi 

u9 

2.92 

0.28 
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8.5 
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MS$4 

0.18 
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1 .4 

31,1*91 
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TABLE  3-3:  The  distribution  function  of  minimal  distance  between 

cloud  pairs  using  a  log-normal  function 


aJ  -  P3  exp  {-  ( InO  -  P1)2/2P22} 


Location  and  center 
coordinates 


Okinawa 

N27°20 ' E129°07 1 


-0.478 
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0.556 


272.748 
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-0.560 
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Pacific  ocean,  NW  of 
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-0.639 


0.529 

0.522 
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370.676 


Ocean  East  of  Florida 
N23°04'W078o44' 


-0.745 

-0.770 


0.526 

0.548 


138.794 
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Atlantic  Ocean  near 
Cape  Verde  Islands 
N15°50iW24°24i 


-0.829 

-0.843 


0.604 

0.590 
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TABLE  3.4:  Examples  of  fit  to  Weibull  function 


■  aOn  exp  {-bDn+1) 


of  the  distance  D 

between 

nearest 

neighbours 

Location  and  center 
coordinates 

n 

b 

Ocean  West  of  Israel 

4 

1.53 

1 .79 

221.50 

N3 10z*6  *  £34°1 1 1 

7 

1.40 

1 .60 

187.69 

Pacific  Ocean  near  Okinawa 

4 

1 .20 

■SSI 

631.36 

N27°20 1 E 1 29°07 1 

7 

1.22 

wxm 

652.90 

Pacific  Ocean  near  Marshall 
Islands 

4 

1.28 

N1 5°50 'W179°20 ‘ 

7 

1.31 

Atlantic  Ocean  near 

Cape  Verde  Islands 

4 

0.83 

1.72 

142.1  16 

N15°50 1 W024°24 1 

7 

0.86 

1.91 

144.543 

Atlantic  ocean  near 

Cape  Verde  Islands 

4 

1.01 

1.95 

217.382 

Nl4o24'W024°45' 

7 

1.21 

2.19 

242.175 

! 
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;  Chapter  4:  On  the  fractal  dimension  of  broken  Cumulus  cloud  fields 

4 . 1  INTRODUCTION 

Clouds  are  important  moderators  of  the  Earth's  radiation  balance  on 
|  both  regional  and  global  scale.  They  are  also  symptoms  or  diagnostic 

features  in  experimental  and  theoretical  studies  of  both  weather  and  climate 
thermohydrodynamics . 

I'any  of  the  remote  observations  of  the  climate  and  weather 
system  necessary  for  the  implementation  or  validation  of  models  or  theories 
are  strongly  influenced  by  the  presence  of  clouds. 

Clouds  are  one  of  the  facets  of  the  problem  of  the  parameterization 
of  sub-resolution  phenomena  for  weather  prediction  and  climate 
modelling.  Finally,  the  morphology  of  cloud  fields  is  important  for 
the  planning  and  analysis  of  all  remote  sensing  data  of  planetary 
atmospheres  and  surfaces  (e.g.  Lovejoy  and  Schertzer,  1986). 

These  are  the  main  reasons  for  our  undertaking  an  extensive  and  intensive 
analysis  of  clouds  in  high  resolution  LANDSAT  satellite  data.  These  will  be- 
in  the  future  -  augmented  by  SPOT,  aircraft  and  other  available  high-resolution 
data  on  cloud  field  morphology. 

Because  of~the  fractal  nature  of  cloud  shapes,  knowledge  of  the  fractal 


dimension  or  dimensions  of  cloud  fields  is  fundamental  to  the  accurate 
description  and  parameterization  of  many  cloud  effects  in  applications  varying 
from  climate  modelling  to  remote  sensing. 
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4.2  DATA  analysis 


Determination  of  the  fractal  dimension  Involves  the  analysis  of  the  area 
and  perimeter  of  each  discernible  cloud. 

The  first  problem  is  then  the  determination  of  a  lower  limit  of  brightness 
of  a  cloud  pixel  over  a  variable  background.  This  difficult  task  is  one  of  the 
reasons  it  was  decided  to  analyse  cloud  data  first  over  the  ocean  where  cloud 
contrast  is  higher  (Joseph  and  Cahalan,  1986,  Joseph  et  al,  1986). 

M;ny  different  techniques  were  tried  based  on  one  and  two-dimensional 
histograms  as  well  as  other  types  of  statistical  analysis  e.g.  spatial  coherence. 
The  technique  finally  adopted  is  described  elsewhere  and  based  on  comparison  of 
size  distributions,  their  moments  and  "noise"  (clouds  less  than  4  pixels  in 
area)  on  a  man-computer  interactive  basis  (Joseph  and  Cahalan,  1986,  Joseph  et  al, 
1986) . 

Secondly,  there  is  the  problem  of  determining  the  length  of  the  perimeter 
of  a  cloud.  This  is  especially  important  for  clouds  with  "holes"  and  even 
more  so  for  the  smaller  clouds. 

We  have  carried  out  a  study  of  three  methods  of  calculating 
the  cloud  perimeter  in  the  course  of  which  one  image  has  been 
analysed  in  great  detail  by  our  U.S.  collaborator. 

The  main  aim  was  a  study  of  the  effect  of  different  choices  of  evaluations 
of  the  perimeter  length  on  the  value  of  the  fractal  dimension.  Three  main 
methods  were  used  to  calculate  the  perimeter  -  counting  perimeter  pixels, 
summing  the  length  of  all  edges  of  pixels  facing  out  from  a  cloud, and,  thirdly, 
method  two  with  an  additional  cutting  of  corners.  All  three  methods  give  fairly 


similar  results  for  the  slope  of  the  log  P,  log  A  curve  for  all  but  the  smallest 
clouds  (less  than  12  pixels  in  area  or  r^<  0.1  km.  In  the  case  of  the  latter,  the 

first  two  again  give  similar  results.  The  last  method  depends  on 
correcting  not  only  the  perimeter  but  the  area  as  well  and  also  on  proper  filtra¬ 
tion  of  the  very  smallest  clouds  -  those  less  than  4  pixels  in  area. 

We  adopted  the  first  method  in  our  study  but  did  not  take  into  account  clouds 

of  area  of  three  pixels  and  less  as  these  are  within  the  noise  level  of  the  data. 

Thirdly,  there  is  the  question  of  taking  into  account  the  fact  that  some 

pixels  may  be  only  partially  cloud-filled.  This  is  particularly  important  for 
the  small  clouds.  It  is  also  connected  to  the  first  problem  mentioned  in  this 
section,  namely  the  identification  of  the  cut  between  cloud  and  background 
intensities.  After  trying  several  approaches  (e.g.  Joseph  et  al,  1986)  we  feel 
that  further  reducing  the  number  of  cloud  pixels  already  identified  by  the  cut, 
by  applying  some  assumption  on  the  perimeter  pixels  -  e.g.  fifty  percent  cloud 
cover  -  is  artificial  and  redundant.  Moreover,  as  shown  in  the  previous 
chapter,  the  derived  properties  are  not  sensitive  to  changes  in  the  location  of 
the  cut  of  the  order  of  one  or  two  grey  levels.  Therefore  any  pixel  identified 
as  a  cloud  pixel  by  our  iterative  analysis  was  counted  as  part  of  the  area  or 
perimeter  of  a  cloud. 


4.3  THE  FRACTAL  DIMENSION  OF  CLOUD  PERIMETERS 


We  have  put  particular  emphasis  in  this  paper  on  the  relationship  between 
areas.  A,  and  parameters,  P,  of  clouds.  A  previous  study  on  such  properties 
of  clouds  fields  (Lovejoy,  1982)  showed- through  the  existence  of  a  constant 

h 

slope  of  (log  P  -  log  A)  curves  in  a  range  of  effective  cloud  radius  (re*  (A/tt) 
between  1  and  1000  km  the  validity  of  a  power  law,  connecting  the  area  and 
the  perimeter  in  the  above  range  of  sizes.  A  second,  very  recent  study  showed 
the  existence  of  two  slope  regimes  (Rys  and  Waldvogel  1986) . 

One  possible  meaning  of  these  results  has  been  discussed  in  the  literature 
(Mandelbrot,  1982)  in  terms  of  the  fractal  dimension  of  the  3-D  cloud  field  and 
of  that  of  a  2-D  cut  through  the  latter  -  e.g.  a  satellite  image  of  a  cloud  field 
the  slopes  of  spectra  of  atmospheric  turbulence  and  scales  of  motions  may  be 
indicated  (ibid  1982). 

The  use  of  fractal  dimensions  to  describe  the  geometry  of  continuous  but  non 
differentiable  curves  and  surfaces  has  been  recently  introduced  (e.g.  Mandelbrot, 
1982),  Clouds  are  such  self-similar  "scaling  fractals"  e.g.  objects  with  no 
intrinsic  scale,  having  a  perimeter  which  depends  only  on  the  resolution  and  the 
"fractal  dimension"  of  the  object.  Empirically,  in  our  case  this  means  that  a 
satellite  image  of  a  cloud  field,  taken  over  a  uniform  dark  ocean  background  with 
no  evident  geographical  features,  would  be  statistically  indistinguishable  from 
an  image  of  any  small  portion  of  the  scene  if  it  had  the  same  number  of  pixels 
as  the  original. 
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Physically*  it  may  mean  that  there  exist  no  preferred  scales  of  motion 
in  the  size  range  where  th.  relationship  applies.  This  may  have  Important 
consequences  for  several  theoretical  and  applied  fields  of  research  in  the 
atmospheric  sciences . 

Lovejoy's  data  are  shown  in  Figure  4.1.  This  is  a  log-log  plot  of  the  area 
versus  the  perimeter  of  rain  and  cloud  areas  determined  by  a  threshold  technique 
from  radar  and  satellite  data.  The  lower  solid  line  is  a  least-squares  fit 
implying  a  power  law  of  the  form 

2  j  j 

area  ***  (perimeter) 

where  d  is  the  "fractal  dimension"  given  of  the  perimeter,  d  =  1.35.  If 
cloud  perimeters  were  smooth  curves  (i.e.  one-dimensional)  then  they  would 
enclose  an  area  proportional  to  the  square  of  the  perimeter,  so  that  the 
slope  would  have  been  2 .  This  case  is  shown  here  for  comparison  as  the  upper 
solid  line  labelled  d  ■  1 . 

Note  that  a  synoptic-scale  system  would  lie  at  the  upper  extreme  of  this 
plot.  A  cloud  covering  one  kilometer  In  diameter,  the  limit  of  resolution  of 
present  meteorological  satellites,  would  be  at  the  lower  extreme.  Over  the 
extensive  sub-tropical  high  pressure  systems  of  the  central  ocean?  and  also 
elsewhere  one  finds  persistent  fields  of  cumulus  clouds  most  of  which  have  a 


diameter  well  below  one  kilometer.  Such  sub-resolution  cloud  fields  (for 
meteorological  satellites)  contribute,  inter  alia,  to  errors  in  satellite- 
based  estimates  of  cloud  cover  and  sea-surface  temperature.  The  idea  that 
clouds  are  scaling  fractals  suggests  that  effects  of  the  smaller  scales  on  are 
area-averaged  quantities  could  be  estimated  from  the  large-scale  data  themselves. 
However,  we  shall  see  that  the  process  of  convection  and  the  existence  of  a 
planetary  boundary  layer  introduce  a  second  fractal  dimension  which  only  becomes 
evident  at  the  higher  resolution  available  from  Landsat. 


Figures  4.2  (a),  (b) ,  (c)  and  (d)  show  a  set  of  (log  P,  log  A)  diagrams  in 
the  blue  spectral  band  (MSS4) ,  The  other  band  analysed  exhibits  identical 
results  and  is  therefore  not  shown  here.  All  results  in  both  bands  are  available 
in  the  Appendices. 

The  cloudiness  of  the  images  analysed  in  the  figures  varies  from  1.4  to 
32  percent,  which  is  about  the  range  of  our  data.  Most  of  the  clouds  are  small, 
much  below  1  km  in  effective  radius.  This  is  normal  for  the  fair  weather  cumulus 
type  of  cloud  field, as  shown  in  a  previous  paper  (Joseph  and  Cahalan  1986).  The 
numbers  and  letters  show  the  frequency  of  occurrence  of  clouds  at  each  experimental 
point.  The  number  of  clouds  involved  is  very  large  and  given  in  Table  4.1  and  in 
the  Appendices  for  all  images  analysed. 

The  line  with  the  smaller  slope  shows  the  linear  least  square  fit  to  all 
data,  which  emphasizes  the  clouds  below  1  km  in  size. 

It  is  obvious  that  for  clouds  larger  than  a  certain  effective  radius,  this 
fit  does  not  give  the  correct  slope  anymore.  Using  the  deviation  of  the  experimental 
points  from  the  fit  to  all  data,  we  determine  statistically  the  most  likely  point 
for  the  change  in  slope.  The  line  with  the  larger  slope  shows  the  best  fit  to  the 
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experimental  data  above  the  latter  point,  it  is  obvious  that  a  statistically 
significant  increase  in  slope-fractal  dimension  occurs.  This  result  has  been 
corroborated  recently  in  the  literature  ( Rys  and  Waldvogel,  1986). 

All  the  results  described  above  are  shown  in  Table  4.1  for  thirteen 
locations,  scattered  over  the  globe.  The  fractal  dimension  of  the  cloud  field 
is  shown,  as  derived  from  experimental  (log  P  -log  rg)  curves  -  like  Figure  4.1 
and  from  the  slope  of  the  size  distributions,  8  (Joseph  and  Cahalan,  1986). 


The  fractal  dimension  for  all  clouds (fourth  column)  and  for  the  large 
clouds  (r  _>  rbreak^  (slxth  column)  is  shown  for  a  variety  of  locations.  The 
average  value  for  all  clouds  taken  together  is  1.29  with  a  standard  deviation 
of  0.3.  Thi  present  results  for  the  clouds  with  an  effective  radius  larger 


than  r 


break’ 


about  0.5  km,  show  a  fractal  dimension  averaging  1.51  ±  0.12. 


The  value  of  1.37,  close  to  Love joy's  value,  was  found  in  one  case  only. 


The  break  radius,  determined  by  an  analysis  of  the  deviation  of  the  residuals 
of  the  fit  from  a  normal  distribution,  is  shown  in  the  fifth  column  in 
km  and  marked  "Break". 


The  location  of  the  break  in  slope  -  or  fractal  dimension  -  is  on  the 
average  at  0 - s_Wj|  pnd  spreads  from  0.178  to  1.0  km. 

In  column  seven,  we  show  the  fractal  dimension  of  the  perimeter  as  calculated 
from  the  slope,  3  of  the  distribution  of  cloud  radii  (Mandelbrot  1984). 


D  =3-1 
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The  slope  is  determined  from  all  clouds  less  than  a  maximal  radius,  r  ,  which 

max 


is  shown  in  column  8.  The  latter  radius  is  defined  as  the  last  point  that  fits 


the  slope  determined  from  the  majority  of  clouds.  The  tail  of  the  distribution 


contains  a  relatively  small  number  of  large  clouds,  which,  though  not  important 


for  the  statistics,  may  contribute  significantly  to  the  total  cloudiness.  This 


maximal  radius  is  of  the  order  of  0.5  km,  or  a  single  cloud  area  of  about 


0.8  km  .  For  more  details  on  the  analysis,  see  the  other  chapters.  The 


ratio  of  the  fractal  dimension  derived  in  the  latter  way,  D  =  3  -  1,  to  that 


from  the  (log  P  -  log  r  )  curves  is  always  larger  than  one.  There  is  no 


apparent  relationship  between  D  and  the  quantity  $  -  1,  which  averages  at 


6-1  -  1.79  ±  0.22. 


4.4  DISCUSSION 


Lovejoy  (1982)  reached  several  conclusions  in  his  pioneering  recent  paper. 


The  first  is  that  clouds  are  self-scaling  fractal  entities  with  a  fractal  dimension 


of  4/3  of  the  perimeter.  Secondly,  that  clouds  are  fractals  down  to  the  lowest  resol 


able  scale  in  the  atmosphere  and  that  therefore  there  are  no  dominant  scales  in 


the  atmosphere f  except  for  the  radius  of  the  Earth.  Thirdly,  that,  because 


cloud  outlines(or  perimeters) and  their  movements  may  be  assumed  to  be  equivalent 


to  those  of  turbulent  isobars,  there  exists  a  fairly  straightforward  physical 


connection  between  homogeneous  Isotropic  turbulence  in  three  dimensions, 


described  by  the  wellknown  5/3  power  law  for  the  energy  scale  spectrum  and  the 


processes  that  produce  clouds  and  rain. 


The  latter  result  was  theoretically  explained  by  an  improvement  in  the 


theory  of  3D  turbulent  diffusion  in  the  atmosphere  (Hentschel  and  Procaccia  1984). 
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The  first  and  serond  of  these  conclusions  may  be  somewhat  doubtful  because 
of  our  results*  The  fractal  dimensions  that  we  find  are  only  rarely  close  to  4/3 
and  are  usually  significantly  smaller  for  the  clouds  less  than  about  0.5  km  in 
effective  radius.  Secondly,  we  find  that  there  occurs  a  distinct  break  in  the 
slope  of  the  (log  P  -  log  A)  relation  at  all  the  geographic  locations  we  have 
investigated.  This  break  is  at  an  effective  cloud  radius  of  about  0.5  km  and 
varies  between  0.2  and  1.0  km  and  is  corroborated  in  the  literature  (Rys  and 
Waldvogel  1986)  . 

Caution  is  indicated  in  reaching  these  conclusions  because  of  the  possible 
effect  of  holes  on  the  perimeter-area  relationships.  This  effect  of  the  holes 
should  be  larger  for  the  smaller  clouds  and  might  thus  lead  to  the  spurious 
occurrence  of  a  break  in  the  slope  and  to  a  smaller  slope  for  the  small  clouds. 

Therefore  an  analysis  of  the  occurrence  and  size  distribution  of  holes  in 
clouds  was  carried  out. 

The  holes  typically  have  an  effective  radius  distribution  with  a  negative  slope 
larger  than  3  and  range  between  a  minimum  effective  radius  of  44  meters  to  about 
450  meters. 

The  negative  slope  of  the  size  distribution  of  the  holes  is  thus  larger  than 
that  of  the  clouds  encircling  them,  leading  to  rapid  decrease  in  the  possible 
influence  of  holes  with  increasing  cloud  size. 

Furthermore,  the  location  of  the  break  is  always  much  over  our  resolution. 

All  the  results  presented  in  this  paper  take  into  account  the  presence  of  holes 

in  the  determination  of  the  perimeter-area  relationships  (See  Chapter  2).  We 

are  therefore  fairly  confident  that  the  break  (or  possibly  a  functional  dependence 

of  the  fractal  dimension  on  r  )  is  in  fact  real  and  should  be  taken  into  account. 

e 


Fair  weather  cumulus  clouds  usually  have  aspect  ratios  of  the  order  of  one 
(e.g.  Plank  1969),  therefore  the  thicknesses  of  the  clouds  average  1  km  and  vary 
between  0.4  and  2  km.  This  is  about  the  depth  of  the  typical  atmospheric 
boundary  layer  that  supports  the  cloud  formation. 

The  effective  cloud  diameter  at  which  the  change  in  the  value  of  D,  the 
slope  of  the  (log  P  -  log  A)  curve,  occurs,  is  at  about  the  same  size  with  the  same 
range. 

It  may  thus  be  that  this  scale  -  that  of  the  thickness  of  the  boundary  layer  - 
is  the  most  important  scale  of  energy  input  into  the  atmosphere  and  that  the  change 
in  slope  is  due  to  a  change  over  from  3D  turbulence  in  the  planetary  boundary  - 
layer  to  2-D  turbulence  in  the  free  atmosphere  (e.g.  Lovejoy,  1986). 

The  differences  between  the  values  of  the  fractal  dimension  D  above  and 
below  the  break  in  each  case  are  statistically  significant.  The  idea  of  clouds 
being  simple  self-similar  scaling  fractals  down  to  the  smallest  size  is  thus  not 
acceptable  without  modification.  Therefore,  for  instance,  it  is  not  directly 
possible  to  use  the  fractal  dimension  of  an  observed  large  scale,  low-resolution 
cloud  field  from  a  meteorological  satellite  in  order  to  correct  for  the  effect 
of  sub-resolution  cloudiness  on  the  determination  of  sea  surface  temperature  or  of 
surface  albedo.  On  the  other  hand,  climatological  data  on  the  variation  of  fractal 
dimensions  with  cloud  size  or  type  may  make  the  above  practical  applications 
possible. 


5  SUMMARY 


We  have  presented  some  of  the  results  of  an  extensive  ongoing  analysis 
of  LANDSAT  imagery  containing  broken  cumulus  cloud  fields  over  the  ocean. 

The  present  results  are  based  on  a  multispectral  analysis  of  \7  locations. 

It  was  shown  that  two  regimes  exist,  each  with  its  own  fractal  dimension. 
All  Cumulus  clouds  less  than  about  500  meters  in  effective  radius  in  a  cloud 
field  seem  to  display  a  fractal  dimension  of  1.29  ±0.03,  those  larger  than 
500  meters  in  radius  -  a  fractal  dimension  of  1.51  ±  0.12. 

These  results  may  have  important  implications  for  both  modelling  and 
remote  sensing  applications. 

If  cloud  fields  in  fact  are  bi-fractal  in  nature,  it  may  mean  that 
between  the  global  and  micro  scales  of  atmospheric  turbulent  thermo¬ 
hydrodynamics,  there  exists  a  scale  of  the  order  of  the  atmospheric  planetary 
boundary  layer  at  which  energy  is  fed  into  the  larger  scales  from  the  planetary 
surface.  This  may  have  important  applications  to  the  modelling  of  climate. 

In  addition,  if  this  result  turns  out  to  be  universally  true,  it  may 
be  possible  to  corr'-'t  low-resolution  remote  sensing  data  for  the  effect  of 
sub-resolution  clou<-.  s,  using  climatologically  determined  relationships 
between  the  fta^di-  nsions  of  large  and  smaller  scale  clouds. 

Future  work  will  enlarge  the  set  of  locations  studied  over  both  ocean 
and  land,  as  well  as  add  the  analysis  of  ultra-high  resolution  data  at 
different  times  of  day. 
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Figure  4.2: 


A  rnn{;e  n f  examples  of  the  relationship  between 
perimeters  and  areas  for  four  locations. 
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Chap tar  5 ;  Application  to  Cumulus  Cloud  Convection 


5.1  The  mean  mass  flux  at  cloud  bottom  altitude 


An  important  application  of  our  present  data  to  atmospheric  physics  is 
an  attempt  to  estimate  the  mean  mass  flux  at  cloud  bottom  altitude.  Ue 
assume  that  the  clouds  in  a  fair  weather  broken  cumulus  field  are  non- 
interacting  transient,  growing  and  decaying,  convective  plumes  with  a 
power- law  size  distribution. 


dN(r) 

dr 


and. cloud  mass  flux  m(r),  both  independent  of  altitude.  Cloud  sizes  are 

bounded  between  a  minimum  effective  radius,  r  .  ,  and  maximum  radius,  r 

min  max 

Under  these  circumstances,  the  mass  flux  per  unit  area  of  the  cloud 
field  at  cloud  bottom  altitude  is  given  by 


dr'  dN(r 1 )  m(r')  /F 


where 


m(r)  -  r 


1x10^  kg  m  ^s  *  disturbed  conditions 


3a  =  0.75 


and  F  is  the  magnitude  of  the  area  investigated.  This  analysis  follows  a  model 
o 


presented  in  the  literature  (Fraedrich,  1976,  1977.') 


Integration  results  in: 


0-2  3a  '■ 


max  o 


The  area  A,  and  accumulated  number  N  of  clouds  can  be  similarly  expressed 


IE  (r3"8  -  r3"8  ) 

3-0  '■max  C  ; 


JL  /r1-8  -  r1'8  ) 

1-0  '•max  ’ 


All  three  accumulated  quantities  can  be  expressed  by  a  constant  and  a 
profile  function  in  the  following  manner. 


M  m  (r) 
o 


N  n  (r) 
o 


A  a  (r) 
o 


where 


Mo  ’  F2  -55"  (rmln  -  rmlx>  /Fo 


m* (r)  -  (r2-8-  r2'8  )/  (r2’6  -  t2 -8) 
max  min  max 


»'■  i  -  ~  •  -.**  .‘*y  k"*  *.*•  / 
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and  M  ,  N  and  A  are  the  values  of  M(r) ,  N(r)  and  A(r)  at  r  =  r  .  . 
o  o  o  min 


The  mean  mass  flux,  Mo>  can  be  expressed  in  terms  of  the  observed  quantities 

N  ,  A  ,  3,  r  .  ,  r  and  F  ,  namely 
o  o’  ’  min’  max  o’  J 


.  ,  2-6  _  2-3 

1  1  (3-1)  (3-3)  min  max 


o  o  (4 


The  units  of  K  are  km 


,  r  .  and  r  are  in  km,  F  in  km 
min  max  o 
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and  Mq  in  kg  m  s  is  then  divided  by  10  to.  get  the  values  shown  in 
Table  5.1,  which  are  given  in  terms  of  increasing  percentage  of  cloudiness 


and  compared  with  three  representative  values  taken  from  the  literature. 


The  results  are  also  shown  in  Fig. 5.1  which  shows  the  mean  mass  flux 
at  cloud  bottom  as  a  function  of  cloudiness  in  percent.  A  least  square 
linear  fit  was  made  between  the  computed  mean  mass  flux  and  the  percentage 
of  cloudiness.  The  correlation  and  fits  are  fairly  good  considering  the 
fact  that  our  data  are  basically  "snapshots"  derived  from  seventeen  different 
locations  at  different  seasons  at  different  mid-morning  times  and  under 
different  synoptic  conditions. 

The  numerical  regression  equation  is 

M0  =  0.59  +  0.049  *  P  _  (5) 

-2  -2  -1 

where  is  in  10  kg  m  s  ,  P  is  in  percent  and  the  correlation  coefficient 
is  0.78. 

The  mean  mass  flux  thus  seems  to  be  linearly  proportional  to  cloudiness 
as  follows  from  our  equations  (3d)  and  (3h)  as  well  as  from  experimental 
data  • 

As  already  mentioned  above,  our  computed  fluxes  represent  an  application 
of  a  particular  transient  plume  theory  of  convection  to  our  size  distribution 
data,  assuming  conditions  to  be  disturbed,  =  1x10  kg  m  s  (Fraedrich 
1976,  1977).  The  size  distribution  data  are  collected  from  many  different 
locations  and  seasons  at  mid-morning  times.  The  slightly  low  correlation  of 
the  fit  (0.78)  is  therefore  to  be  expected.  At  all  locations  the  presence  of 
clouds  probably  represents  fair  to  disturbed  weather  conditions. 
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The  GATE  A/B  results (Houze  and  Cheng  1981)  also  represent  convective 
conditions,  but  averaged  over  the  whole  GATE  area  for  A/B  conditions  which 
include  Fair  Weather  Cumulus  as  well  as  more  highly  convective  situations. 

The  fact  that  the  GATE  A/B  area  fluxes  are  slightly  higher  than  ours  for 
the  same  average  amount  of  cloudiness  is  therefore  understandable. 

Fraedrich's  data  for  50%  cloudiness  are  slightly  lower  than  our  prediction  for 
such  cloudiness.  Again  this  is  understandable  as  they  represent  ATEX 
undisturbed  conditions  only. 

Although  the  agreement  of  the  fluxes  computed  from  our  data  appear 
fairly  well  to  the  in-situ  results,  our  results  should  be  only  taken  as  an 
indication  that  further  studies  are  strongly  indicated. 

It  should  be  mentioned  in  this  context  that  one  of  the  conclusions  from 
the  GATE  experiment  is  the  importance  of  the  few  large  cloud  areas  for  rain¬ 
fall  and  cloudiness.  This  may  indicate  that  one  needs  to  deal  with  multi¬ 
modal  size  distributions  . 

If  one  assumes  that  the  mass  fluxes  of  individual  non-interacting  clouds 
do  not  vary  with  height,  one  can  calculate  the  change  with  height  of  the 
average  mass  flux  of  the  bottom  of  the  clouds  directly  from  the  profile 
function  (3e).  type  of  clouds  we  have  analyzed  have  aspect  ratios  - 

diameter  to  thickness  ratios  -  of  the  order  of  one  (Plank  1969).  Therefore, 
we  can  transform  the  dependence  on  r  to  a  dependence  on  height  h  above  the 
common  cloud  bottom, by  using  that  fact,  e.g. 


2r 

h 


1 


(*) 


The  mass  flux  then  decreases  approximately  as  the  ratio  of  the  total  cloud 


area  at  height  h  ■  2r  to  that  at  cloud  bottom. 


r(3-6) 

max 

max 


-  (h/2) 


(3-0) 

min 


(3-6) 


(7) 


The  function  peaks  at  the  bottom  of  the  cloud  layer,  and  its  decrease  to 
zero  above  that  height  depends  on  B  and  r 

max 

The  larger  r  is,  the  slower  will  the  profile  function  decrease  with 
height.  The  dependence  on  6  is  weaker,  as  the  latter  varies  in  our  data  between 
3-3  and  3.5. The  importance  of  taking  into  account  all  clouds  is  thus  clear. 
However,  our  size  distribution  combined  with  the  sample  parameterization  model 
for  the  flux  cannot  reproduce  the  maximum  in  the  profile  at  about  450  mb  in 
the  GATE  A/B  data,  which  is  probably  due  to  meso-scale  phenomena,  not  taken 
into  account  in  the  model  used  by  us  (Houze  and  Cheng  1981) .  The  rapid  decay 
of  MQ(h)  with  height  above  cloud  base  shown  by  Fraedrich  (Fraedrich  1977)  is 
reproduced  approximately. 

5.2  The  ratio  of  updraft  to  downdraft  in  convective  Cu  clouds 

Another  application  of  our  data  to  the  convection  problem  is  an  estimate 
of  the  ratio  of  average  updraft  speeds  to  average  downdraft  speeds.  If  one 
assumes  that  each  active  cumulus  cloud  is  isolated  from  all  others  and  that 
mass  is  conserved,  then,  from  mass  continuity, we  must  have  a  definite  ratio 
between  the  area  of  the  updraft  and  that  of  the  downdraft  for  a  given  cloud.  We 
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tried  to  approximate  that  ratio  in  the  following  way,  assuming  the  total 
area  of  each  visible  cloud  to  represent  an  updraft. 

We  calculated  for  each  cloud  of  given  radius  r,  the  nearest  distance,  D, 
from  those  neighbours  whose  areas  are  at  least  one  tenth  greater.  The  latter 

condition  was  decided  upon  in  order  not  to  calculate  trivial  distances 
between  a  cloud  and  a  small  puff  close  by  that  may  be  part  of  it.  We  then 
plot  the  ratio  r  /(CD/2)  -  r  )  which  we  term  y.  The  latter  ratio  should  be  an 

approximation  to  the  ratio  of  cloud  area  to  clear  area  around  it  that  is  part 
of  the  same  convection  cell.  Only  one  half  of  the  nearest  distance  was  taken 
for  the  same  reason. 


Typical  results  are  shown  in  Figures  5.2(a),  (b)  and  (c) .  These  are 

two-dimensional  plots  of  the  ratio  y  versus  the  effective  radius  r  .  Because 

e 

of  the  very  large  number  of  clouds  involved  in  each  case  -  of  the  order  of 
thousands  -  the  points  run  together  and  form  a  pattern,  similar  for  all  cases 
analysed.  This  pattern  has  the  following  properties.  First,  the  great 
majority  of  points  cluster  in  the  upper  right  (positive)  quadrant  of  the 
(y,  rg)  plane  in  the  vicinity  of  the  origin.  There  are  no  negative  values 
of  y,  except  the  rare  cases  when  two  elongated  clouds  are  very  close  together. 
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This  result  indicates  that  the  mode  or  preferred  value  of  y  is  very 
small  and  positive.  Secondly,  on  analysis  of  the  figures,  it  becomes  clear 
that  this  mode  value  is  independent  of  the  effective  cloud  radius. 

In  order  to  study  these  results  in  some  more  detail,  we  show  in 
Tables  5.2(a)  and  5.2(b)  the  same  data  for  both  bands  (MSS4  and  MSS7)  of 
Case  7072. 

Because  of  the  large  number  of  clouds  involved,  1 1 , 300 , sp anning  a  range  of 

2 

areas,  A,  from  0.002  to  20  km  ,  the  effective  radii 


re  =  (A/it)  4  (8) 

where  divided  into  groups  0f  equal  size  intervals  when  there  are  many  clouds 

and  the  last  10  clouds  in  one  final  interval.  In  this  way  the  few  large  clouds  can 

still  show  their  own  behaviour  w.r.t,  those  smaller. 

The  number  of  cases  in  each  radius  range  is  thus  not  equal.  However,  due 
to  the  large  number  of  clouds  involved  in  each  range,  except  for  the  largest 
clouds,  this  is  not  a  serious  problem.  We  find  that  in  the  cases  considered, 
the  ratio 

2  2  9 

r  /((D/2)  -  r  ) 

again  shows  the  two  above  properties,  which  are  typical  of  all  cases  analysed. 


First,  the  peak  value  of  the  histogram  of  the  ratio  in  each  size  range 
is  between  0  and  0,025.  The  majority  of  ratios  are  in  this  range  as 
shown  in  Tables  5. 2. a  and  5.2.b. 

Secondly,  the  peak  value  is  almost  independent  of  cloud  size  up  to  and 
including  the  largest  clouds. 

It  seems  to  us  that  the  properties  indicated  by  this  result  are  that  the 
clouds  are  formed  in  isolated  cells  in  a  region  of  fairly  homogeneous 
horizontal  distribution  of  water  vapor  content.  Turbulent  convergence  of 
water  vapor  flux  led  to  the  formation  of  a  cloud  in  each  cell  concentrating 
the  water  of  the  latter  into  the  visible  cloud  volume.  The  area  or  volume  of  each 
turbulent  cell  determines  the  amount  of  available  water  vapor  and  thus  the 
size  of  the  cloud.  The  ratio  of  water  vapor  source  region  to  cloud  region  for 
each  cloud  should  be  constant  and  lead  to  our  result. 

The  range  of  values  found  by  us  may  be  compared  to  a  theoretical  estimate 
of  the  ratio  of  a  cloud  radius  to  its  "radius  of  influence"  .  The  above 
ratio  may  be  approximated  by  r(D/2)  and  using  our  data  gives  a  range  of 

0  <  r/ (D/2)  <  0.16 

as  compared  to  the  theoretical  value  of  about  0.2  for  amounts  of  cloudiness 
between  5  and  15  percent  (Soong  and  Ogura  1973,  Cho  1978). 
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Figure  5.1:  ~  The  dependence  of  the  mean  mass  flux  at  cloud 

bottom  as  a  function  of  cloudiness  in  percent 
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Figure  5.2: 


T%e  parameter  V  as  a  function  of  r  for  three  cases: 

(a)  Case  7067,  cloudiness  1.47;  (b^  Case  7025,  cloudiness 
13*5  .  (c)  Case  7024  .  cloudiness  33%.  All  cases  in  Band  MSS4 

Results  for  MSS7  are  identical. 
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Table  5.1:  Mean  mass  flux  at  cloud  bottom  altitude  in  10* 


-2  -1 


Image  # 

* 

Location 

-2  -2-1 
Mq(10  kgm  s  1 

Cloudiness  Z 

F  km2 
o 

7036 

Ocean  near  Cape  Verde 

Island 

0.113 

0.2 

3.2  (4) 

7034 

Ocean  near  Cape  Verde 

Islands 

0.236 

0.6 

3.2  (4) 

7072 

Gulf  of  Mexico 

0.130 

0.8 

1.8  (4) 

7035 

Ocean  near  Cape  Verde  Island 

0.489 

0.9 

1.4  (4) 

7067 

Gulf  of  Mexico 

0.656 

1.5 

3.2  (4) 

7066 

Gulf  of  Mexico 

1.132 

5.1 

1.9  (4) 

7032 

Ocean  near  Marshall  Island 

1.362 

6.7 

3.1  (4) 

7013 

Med.  W.  of  Israel 

1.796 

7.5 

6.2  (3) 

7029 

Ocean  E.  of  Florida 

1.329 

8.7 

1.6  .<4) 

7064 

Ocean  E.  of  Florida 

0.522 

12.1 

2.0  (4) 

7025 

Ocean  E.  of  Florida 

1.933 

13.2 

1.7  (4) 

7033 

Pacific  Ocean  in 

Tuamotu  Archipelago 

0.812 

14.3 

3.2  (4) 

7028 

Ocean  E.  of  Florida 

1.344 

17.4 

1.7  (4) 

7026 

Ocean  E.  of  Florida 

1.884 

18.8 

1.5  (4) 

7030/31 

Okinawa 

2.181 

27.0 

2.2  (4) 

7024 

Ocean  E.  of  Florida 

2.083 

33.1 

8.0  (3) 

7060 

Bahamas 

1.967 

37.4 

l.9  (4) 

Cheng  (1980) 

1.670 

10.0 

GATE  A/B  average 

Houze  and  Betts  (1981) 

2.000 

15.0 

Ocean  E.  of  Florida 
(Fraedrich  1977) 

3.000 

50.0 

it 

For  exact 

coordinates,  path  and  row  -  see 

Appendix  I 

V-<, 


-  127  - 


TABLE  5.2a;  2-D  Histogram  of  the  ratio  of  cloud  area  to  cloud  area 

—  of  Influence  as  a  function  of  cloud  effective  radius  Y  =  — 
Location:  Gulf  of  Mexico,  ID  #7072,  MSS4  ( 


2  2 
(D/2)  V 
e 


0.575-0.600 

2 

2 

1 

0 

0 

0 

0.550-0.575 

1 

1 

0 

1 

0 

0 

0.525-0.550 

1 

1 

0 

1 

0 

1 

10.500-0. 525 

0 

2 

0 

1 

0 

0 

b.  475-0. 500 

3 

3 

0 

0 

0 

0 

0.450-0.475 

1 

1 

0 

0 

0 

0 

0.425-0.450 

1 

0 

3 

1 

0 

0 

0.400-0.425 

1 

2 

3 

1 

0 

0 

0.375-0.400 

1 

0 

1 

0 

1 

0 

0.350-0.375 

2 

1 

0 

0 

1 

0 

0.325-0.350 

2 

4 

1 

1 

0 

0 

0.300-0.325 

2 

1 

0 

0 

0 

0 

0.275-0.300 

3 

1 

0 

0 

1 

1 

0.250-0.275 

3 

5 

1 

1 

0 

0 

0.225-0.250 

4 

6 

1 

2 

1 

0 

■0.200-0.225 

7 

8 

3 

0 

1 

0 

3.175-0.200 

7 

3 

2 

0 

1 

0 

3.150-0.175 

10 

7 

4 

0 

0 

1 

p. 125-0.150 

7 

16 

5 

0 

2 

1 

!0. 100-0. 125 

7 

12 

3 

5 

0 

2 

0.075-0.100 

14 

17 

3 

i 

1 

0 

0.050-0.075 

24 

25 

7 

i 

2 

1 

0.025-0.050 

51 

-n 

8 

0 

0 

1 

0.000-0.025 

54 

37 

10 

3 

0 

0 

Introduction  to  the  Appendices  I  -  IV 


Appendices  I,  II,  III  and  IV  summarize  in  concise  form  a  great  part 
of  the  basic  information  acquired  during  the  course  of  this  study. 

Appendix  I  identifies  the  analysed  images  in  several  ways.  In 
addition,  information  is  provided  on  the  parts  of  each  image  that  were 
not  analysed  and  the  reasons  for  this,  which  are  of  interest  in  themselves. 
There  are  mainly  three  . 

The  first.  A,  is  the  presence  in  the  image  of  land  surfaces  or  shallow 
ocean.  It  is  impossible  to  distinguish  reliably  between  clouds  and  such 
variable  backgrounds  even  with  intensive  man-computer  interaction.  Con¬ 
sequently,  we  have  not  analysed  those  parts  of  the  images  mentioned  in  the 
appendix.  In  addition,  we  will  refrain  from  analysing  any  images  over  land, 
until  appropriate  techniques  can  be  developed,  which  will  probably  have  to 
include  utilisation  of  information  on  the  surface  reflectivity  and  its 
spatial  distribution. 

The  second  reason,  B,  is  the  presence  of  fog  or  haze  over  the  image 
with  an  intensity  close  to  that  of  the  smaller  clouds.  This  makes  it  wellnigh 
impossible  to-discriminate  a  great  number  of  the  smaller  clouds,  part  of  the 
studied  population,  against  the  background. 

The  last  reason,  C,  is  the  partial  absence  of  an  image  or  errors  in 
the  grey  levels  in  the  data  base.  It  may  also  happen  that  there  are  over¬ 
estimates  of  cloudiness,  especially  in  cases  of  reported  cloudiness  of 
ten  percent  or  more  in  the  acquired  satellite  data. 


The  second  Appendix,  II,  summarises  for  the  two  bands  analysed,  4  and  7, 
the  area  analysed  in  pixels,  the  criterion  used  for  the  grey  levels,  in  each 

case,  to  distinguish  between  cloud  and  background,  the  exact  percentage  of 

cloudiness  of  the  area  analysed.  Moreover,  we  specify  the  fractional  area 
of  clouds  with  areas  of  less  than  4  pixels  as  "noise",  as  well  as  the  number 
of  such  small  clouds  versus  the  number  of  clouds  of  larger  size. 

Appendix  III  summarises  all  available  empirical  fits  to  cloud  size 

distributions  in  the  form  of  power  laws.  Both  the  fit  up  to  a  maximal 

effective  radius  as  well  as  that  using  all  clouds,  are  given,  showing  the 
probable  need  for  bi-modal  distributions. 

Secondly,  all  (log  P,  log  A)  fits  given  again  for  each  histogram  up  to 
a  maximal  size  as  well  as  using  all  clouds.  The  radius  of  the  break  is 
identified  in  each  case. 

In  the  last  column,  the  two  average  fractal  dimensionsof  the  perimeters 
in  each  case  are  given. 

Appendix  IV  summarises  fits  to  the  distribution  functions  of  the  minimal 

distance  between  cloud  pairs.  Both  log-normal  and  Weibull  distributions  are 

2 

used.  All  distributions  are  normalised  to  100  km  ,  to  make  comparisons 
possible  between  the  various  cases. 

Some  additional  general  information  is  required  for  using  the  appendices. 

p 

The  area  of  one  pixel  is  0.003249  km.  "Large  clouds"  are  those  clouds  that 
have  an  area  of  more  than  three  pixels.  Only  these  are  considered  clouds.  "Small 
clouds"  are  those  of  less  than  4  pixels  in  area.  This  is  the  limit  of  our  resolu¬ 
tion  and  consequently  the  "small  clouds"  are  considered  to  be  noise: 

There  are  2983  lines  and  3548  columns  in  a  full  image. 
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Areas  of  Image  Not  Processed 
R  *  Reason  for  Not  Processing: 

Image  Image  Place  Path  Date  A  =  Land  or  Shallow  Ocean 

No.  ID  Row  YY  DOQ  B  =  Fog  Over  Ocean 

C  *  Data  Absent  or  in  Error 


Lines  Columns  R  Lines  Columns  R 


7013 

2555-07291 

Israel 

188-038 

76-212 

7024 

22072-15003 

Bahamas 

014-044 

80-268 

7025 

22073-15055 

Straits 

FI orida 

015-043 

80-269 

7026 

22056-15113 

Straits 

Florida 

016-043 

80-252 

7028 

22109-15054 

Straits 

Florida 

015-043 

80-305 

7029 

22127-15052 

Straits 

Florida 

015-043 

80-323 

7030- 

7031 

30277-01144 

Okinawa 

78-341 

7032 

30346-21395 

Paci fic 
Ocean 
near 
Marshal  1 
Islands 

083-049 

79-045 

7033 

30369-18510 

Paci fic 
Ocean 
Toamotu 
Islands 

052-071 

79-068 

7034 

30364-11203 

Cape 

Verde 

W. Africa 

226-050 

79-063 

7035 

30364-11201 

Cape 

Verde 

W. Africa 

226-049 

79-063 

7036 

30364-11194 

Cape 

Verde 

W.  Africa 

226-0-18 

79-063 

7060 

22036-14590 

Bafcems 

014-041 

80-232 

7064 

22038-15094 

Ocean  E 
of 

Florida 

016-039 

80-234 

7066 

22312-15310 

Gulf 

Mexico 

020-040 

81-1-12 

7067 

30661-15303 

Gulf 

Mexico 

020-040 

79-360 

7072 

31565-15595 

Gulf 

Mexico 

024-040 

80-168 

1-  700 

900-3548 

A 

701-1200 

1200-3548 

A 

1201-1500 

1600-3548 

A 

1501-2100 

1900-3548 

A 

2101-2983 

1-3548 

A 

1-1200 

1-950 

A 

1-2983 

2650-3548 

A 

2100-2983 

750-2200 

A 

1-1800 

1-1200 

1-1400 

1201-3548 

1-1200 

1-  950 

A 

1-2983 

2650-3548 

A 

2100-2983 

750-2200 

A 

1-1200 

1-  700 

A 

1-1500 

2400-3548 

A 

1501-2983 

2600-3548 

A 

2100-2983 

500-2100 

A 

1-  219 

1245-1700 

A 

407-  749 

934-1330 

A 

1075-1277 

544-  950 

A 

1970-2340 

1-  650 

A 

No  holes 

1100-2250 

1520-2200 

A 

2400-2679 

200-  700 

A 

880-  930 

1-3548 

C 

2680-2800 

1-3548 

C 

450-  700 

1500-1900 

A 

100-  700 

1900-2500 

A 

1-1299 

1-3548 

B 

1300-2983 

2601-3548 

B 

2200-2983 

1-1320 

A 

2640-2983 

1320-2200 

A 

1980-2200 

1-  200 

A 

’900-2983 

1-2200 

A 

2420-2983 

2201-3548 

A 

660-1500 

2200-3548 

B 

1760-2983 

1400-3548 

B 

1  - 1800 

1-3548 

B 

No  holes 

1-2983 

1-1030 

C 

1-  220 

1320-1685 

B 

1-  880 

1686-3548 

A 
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Image 

No. 

MSS 

Area 
Analysed 
in  Pixels 

Grey  Level 
Threshold 
for  Clouds 

Percentage 

of 

Cloudiness 

Percentage 

of 

Noise 

No.  of 

Large 

Clouds 

No.  of 

Small 

Clouds 

7013 

4 

190  9440 

20 

8.5922 

0.4038 

1269 

- 

7 

190  9440 

3 

7.3577 

0.3482 

1242 

7024 

4 

2473345 

21 

32.9202 

0.1037 

2521 

2566 

7 

2473345 

3 

33.1866 

0.1156 

2420 

1728 

7025 

4 

5285332 

69 

13.4944 

0.1039 

7329 

3173 

7 

5285332 

2 

12.8875 

0.1315 

7436 

4538 

7026 

4 

4678420 

21 

18.9997 

0.0826 

4738 

2235 

7 

4678420 

3 

18.4689 

0.0793 

4671 

2152 

7028 

4 

5304581 

69 

17.4672 

0.0647 

3349 

2017 

7 

5304581 

2 

17.3171 

0.0742 

3458 

2345 

7029 

4 

4999329 

65 

3.7629 

0.0759 

4235 

2176 

7 

4999329 

2 

8.5214 

0.0930 

4241 

2836 

7030 

4 

6909474 

43 

27.5781 

0.1431 

6836 

- 

7031 

7 

6909474 

3 

26.4643 

0.2021 

7013 

7032 

4 

9650303 

82 

6.6937 

0.0540 

7327 

2976 

7 

9650303 

7 

6.7086 

0.0952 

7548 

6255 

7033 

4 

8276192 

77 

14.5296 

0.0320 

3616 

1794 

7 

8276192 

6 

13.9987 

0.0742 

3562 

4750 

7034 

4 

9716240 

77 

0.5882 

0.0147 

1477 

809 

7 

9716240 

6 

0.6015 

0.0282 

1380 

1634 

7035 

4 

4210868 

81 

0.9156 

0.0311 

1274 

728 

7 

4210868 

7 

0.9053 

0.0403 

1203 

1104 

7036 

4 

9714234 

77 

0.1843 

0 . 0086 

678 

447 

7 

9714234 

6 

0.1981 

0.0749 

709 

5019 

7060 

4 

6685609 

22 

37.6347 

0.1082 

6676 

4146 

7 

6685609 

3 

37.1414 

0. 1402 

6631 

5827 

7064 

4 

6216943 

27 

12.0736 

0.0493 

2440 

1829 

7 

6216943 

5 

12.1883 

0.0485 

2252 

1849 

7066 

4 

5789598 

24 

5.2203 

0.3688 

5172 

13842 

7 

5789598 

5 

4.8984 

0.4109 

5355 

15711 

7067 

4 

9692370 

17 

1.3787 

0.0873 

4450 

5452 

7 
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Fit  to.  Cloud  Radius  Distribution 


Fit  to  (logP,  logA)  Curve 
log10P  =  a  +  blog1QA 
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- 
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all  -0.01273 
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all 
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up  to  r  =1.00  km  -0.34390 
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- 
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all  e  -0.02089 
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all 
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7 
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-2.68639 

- 

7033 

4 

all  -0.65950 
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all 

0.88669 

0.63058 
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up  to  r  =0.56  km  -0.77968 

-2.85170 
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km 

0.89297 
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1.40530 

7 

all  e  -0.55310 
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all 

7 

0.90364 

0.63772 
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-2.82177 
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7034 

4 

all  -0.56379 
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all 

7 
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4 
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from  A=0.50 

km*- 

0.95080 
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7 

all  e  -0.54762 

-1.94276 

all 

7 

0.91381 
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7 

up  to  rg=0.36  km  -1.38271 

-2.75142 
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km 

0.95920 

0.77243 
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4 

all  -0.65404 
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7 

0.83445 
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km 
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0.73822 

1.41267 

7 

all  e  -0.63874 

-1.99064 

all 

7 

0.89696 

0.62973 

7 

up  to  r^‘0.35  km  -1.34412 

-2.98635 

from  A- 0. 50 

km 

0.94904 

0.67445 

7036 

4 
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-0.70365 

-2.82239 

7 

up 

to  r  =0.3  km 

e 

-0.62811 

-2.76809 

7072 

a 

all 

-0.27861 

-1.06982 

4 

up 

to  r  =0.25  km 

-1.68625 

-2.91277 

7 

all  e 
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FIT  TO  MINIMAL  DISTANCE  BETWEEN  CLOUD  PAIRS,  D 


N(D)  -  P3*  exp  -((In  D-P1)2/(2P2)2  ) 

+  *  P3  Computed  for  Normalized  Areas 

N(D)  *  P3*DPl  *exp(-P2*D  ) 

+  =  P3  Computed  for  Normalized  Areas 

Image 

No. 

MSS 

PI 

P2 

P3 

PI 

P2 

P3 
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4 

-0.47792 
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_ 

7 
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-0.77908 

0.55998 

355.6379 

1.20922 

2.18450 

1323.11584+ 
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0.57941 

1 38 . 6066 
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7067 

4 
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2  39 . 07 1 1 + 
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DESCRIPTION  OF  COOPERATION  AND  ASSESSMENT  OF  THE  CONTRIBUTION  OF  BOTH  PARTNERS 


1.  Introduction 


The  work  of  the  two  partners  is  very  tightly  Interconnected  and  complementary. 
In  addition,  however,  each  of  the  two  partners  also  has  separate  aims  and  plans 
w  .r . t .  the  resul ts . 


2.  Contribution  of  U.S.  partner 
The  U.S.  partner  has  provided 

.  catalogues  of  LANDSAT  MSS  data  at  EROS  Data  Center  and  at  NASA-GSFC 

.  60°  swath  of  LANDSAT  4  data  in  the  Pacific  Ocean  is  being  analysed  by  him 

.  two  7-band  scenes  of  LANDSAT  4  TM  data  including  tape  dump  in  the 

Carolinas  area 

.  software  for  reading  LANDSAT  4  TM  tapes 

.  documentation  on  facilities  and  software  available  at  the  Goddard  Space 
Flight  Center 

.  advice  on  Weibull  functions  and  their  use  in  statistics 

.  comparison  of  different  techniques  for  calculating  the  length  of  perimeters 
for  small  clouds 


3.  Contribution  of  Israeli  partner 

The  Israeli  partner  has  provided 

.  size, nearest  neighbour  spacing  and  other  distribution  functions  and  other 
data 

.  non-linear "fits  of  the  above  to  various  probability  distribution  functions 
indicated  by  physical  intuition 

.  application  to  problems  of  convection  in  the  atmosphere 

.  development  of  radiative  transfer  technique  for  cloud  top  height 
determination 

.  software  for  the  analysis,  identification  and  labelling  of  clouds  and 
cloud  fields 

.  choice  of  LANDSAT  3,  4  and  5  imagery  for  analysis  and  their  organization 
into  catalogues 


computer  memory  disk  rental  and  data  purchase  through  independent  funding 


In  addition  to  the  above  contributions,  each  partner  has  pursued  his  own 
specific  plans  with  the  results. 

The  US  partner  emphasizes  analysis  of  the  relationship  of  sub-grid  scale 
convective  processes  on  horizontally  averaged  quantities  in  the  atmosphere  and 
the  large-scale  conditions  in  which  cumulus  convection  occurs. 

The  Israeli  partner's  aims  are  to  use  the  data  in  the  parameterisation  of 
vertical  radiative  transfer  of  energy  in  fields  of  broken  cloudiness  and 
diagnosis  of  conditions  of  atmospheric  convection  dynamics  from  the  cloud  field 
structure,  cloud  height,  etc. 
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Inventory  and  Status  of  Data  -  May  1986 
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